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Abstract 

Using  a  device  to  act  as  a  surrogate  for  atmospheric  turbulence  in  a  labora¬ 
tory  is  necessary  to  build  and  test  optical  systems  for  imaging,  lidar,  laser  weapons, 
and  laser  communications.  Liquid-crystal  spatial  light  modulators  (LC  SLMs)  and 
segmented  micro-electro-mechanical-system  (MEMS)  deformable  mirrors  (DMs)  are 
commonly  used  devices  for  altering  wavefronts  in  order  to  simulate  a  portion  of  at¬ 
mospheric  turbulence.  The  best  location  of  these  devices  was  theoretically  analyzed 
to  obtain  the  broadest  possible  range  of  atmospheric  conditions.  It  was  found  that 
two  phase  screens  should  be  placed  at  the  beginning  of  the  optical  path  to  achieve 
maximum  turbulence  strength  for  an  incident  plane  wave.  It  was  also  revealed  that 
a  layered  model  for  atmospheric  turbulence  strength  can  be  represented  by  the  at¬ 
mospheric  coherence  diameter  that  a  single  segmented  wavefront  control  device  can 
produce.  The  limitations  of  pixelation  effects  on  a  segmented  wavefront  control  device 
were  investigated  theoretically.  The  results  of  this  analysis  were  then  confirmed  in 
simulation.  It  was  found  that  while  LC  SLMs  with  high  bandwidth  have  almost  no 
adverse  affects  from  pixelation,  segmented  MEMS  DMs  have  limitations  as  a  result 
of  the  number  of  mirror  segments  on  a  DM.  The  performance  capabilities  of  several 
available  commercial  devices  are  better  understood  as  a  result  of  this  research. 
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Limitations  of  Segmented  Wavefront  Control  Devices 
in  Emulating  Optical  Turbulence 

I.  Introduction 

This  research  aims  to  find  the  maximum  equivalent  atmospheric  turbulence 
strength  that  can  be  achieved  in  a  laboratory  environment  with  two  phase  screens.  It  is 
revealed  that  a  layered  model  for  atmospheric  turbulence  strength  can  be  represented 
by  the  atmospheric  coherence  diameter  that  a  single  segmented  wavefront  control 
device  can  produce.  The  limitations  of  pixelation  effects  on  a  segmented  wavefront 
control  device  are  investigated  theoretically  and  in  simulation. 

1 . 1  Motivation 

Using  a  device  to  act  as  a  surrogate  for  atmospheric  turbulence  in  a  laboratory 
is  necessary  to  build  and  test  optical  systems  for  imaging,  lidar,  laser  weapons,  and 
laser  communications.  The  United  States  Air  Force  has  several  multi-million  dollar 
programs  that  require  optical  propagation  through  the  atmosphere.  These  include 
the  Airborne  Laser  program,  the  Advanced  Tactical  Laser  program,  and  the  satellite 
imaging  program  at  the  Starfire  Optical  Range  at  Kirtland  Air  Force  Base,  New 
Mexico.  Each  of  these  programs  require  the  propagation  of  light  through  a  turbulent 
atmosphere.  How  atmospheric  turbulence  affects  optical  propagation  must  be  fully 
understood  before  any  of  these  programs  can  try  to  complete  their  respective  missions. 
This  is  because  the  controllers  need  to  fully  understand  those  effects  so  they  can  be 
mitigated  before  any  laser  is  fired. 

Atmospheric  turbulence  effects  are  mitigated  with  the  use  of  adaptive  optics  (AO). 
Large-scale  deformable  mirrors  (DMs)  are  usually  employed  in  these  systems  as  real¬ 
time  AO  devices  that  correct  for  the  effects  of  the  atmosphere.  A  turbulent  wavefront 
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is  one  in  which  the  atmosphere  has  altered  the  phase  of  propagating  light  via  ran¬ 
domly  evolving  distortions  of  the  refractive  index.  Complex  wavefront  sensors  in  AO 
systems  measure  the  amount  of  distorted  phase  added  to  the  light  by  turbulence. 
Then  DMs  are  used  to  correct  for  that  turbulence.  DMs  compensate  by  distorting 
the  shape  of  the  mirror  so  that  the  conjugate  of  the  turbulent  phase  is  applied  to  the 
incident  wavefront.  When  the  turbulent  wavefront  reflects  off  the  DM,  this  conjugate 
phase  flattens  the  wavefront.  There  are  many  areas  of  research  that  enable  this  com¬ 
plex  AO  system  to  operate  correctly  at  real-time  speeds.  Wavefront  sensing,  control 
systems,  and  wavefront  control  devices  are  extensively  researched  in  order  to  improve 
this  process. 

Before  a  DAI  is  put  to  work  in  a  complicated  system  with  real  turbulence,  its  in¬ 
dividual  capabilities  must  be  understood.  These  capabilities  are  first  investigated  in  a 
controlled  laboratory  environment.  Common  laboratory  conditions  are  not  conducive 
to  testing  atmospheric  turbulence.  Therefore,  a  suitable  turbulence  surrogate  must  be 
employed.  Heating  elements  can  be  used  to  create  turbulence  in  a  laboratory,  but  the 
turbulence  they  create  cannot  easily  be  controlled  or  reproduced.  A  better  method  of 
reproducing  atmospheric  turbulence  in  a  laboratory  is  necessary  to  test  the  abilities  of 
wavefront  control  devices.  This  research  focuses  on  the  capabilities  and  limitations  of 
wavefront  control  devices.  Liquid-crystal  spatial  light  modulators  (LC  SLMs),  micro 
electro-mechanical  systems  (MEMS)  DMs,  and  static  atmospheric  phase  wheels  are 
small  wavefront  control  devices  that  are  commonly  used  in  laboratories  for  recreating 
and  correcting  optical  aberrations.  Without  first  testing  a  device’s  abilities,  a  pro¬ 
gram  is  not  able  to  use  that  device  to  its  full  potential.  Failure  of  a  device  to  perform 
in  a  complex  program  could  result  in  failure  of  the  mission  and  significant  monetary 
loss. 

The  solution  to  the  problem  of  emulating  turbulence  in  a  laboratory  lies  in  using 
wavefront  control  devices  to  both  produce  and  correct  for  atmospheric  turbulence.  If 
a  device  is  capable  of  applying  the  conjugate  of  the  turbulent  phase  to  an  incident 
light  source,  then  it  is  also  capable  of  applying  the  turbulent  phase  to  a  plane  wave 
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light  source.  This  means  that  a  simple  coherent  laser  source  can  be  altered  by  a  wave- 
front  control  device  to  recreate  atmospheric  turbulence.  Using  a  wavefront  control 
device  results  in  repeatable  and  controllable  atmospheric  turbulence  parameters  for 
the  distorted  wavefront. 

Another  important  way  to  inexpensively  study  how  the  atmosphere  affects  op¬ 
tical  propagation  in  a  program  is  to  use  computer  simulation.  Wavefront  control 
devices  in  a  laboratory  are  used  to  verify  that  these  computer  simulations  are  accu¬ 
rate.  Without  verification  through  experimentation,  a  simulation  cannot  be  trusted 
to  act  as  a  surrogate  for  atmospheric  turbulence.  Thus,  these  two  techniques  are 
complementary. 

1.2  Goals 

The  primary  objective  of  this  research  is  to  identify  the  best  locations  for  wave- 
front  control  devices  in  a  laboratory  experiment  in  order  to  maximize  the  overall 
turbulence  they  can  produce.  This  research  focuses  on  the  locations  for  two  wave- 
front  control  devices  and  a  plane  wave  light  source. 

To  satisfy  this  goal,  this  research  investigates  the  limitations  in  emulating  atmo¬ 
spheric  turbulence  of  segmented  wavefront  control  devices  such  as  segmented  MEMS 
DMs  and  LC  SLMs.  Another  goal  of  this  research  is  to  identify  and  demonstrate 
the  limits  of  these  devices.  This  goal  is  addressed  through  an  investigation  of  the 
maximum  turbulence  strength  a  segmented  wavefront  control  device  can  emulate 
through  both  theoretical  analysis  and  computer  simulation.  The  maximum  turbu¬ 
lence  strength  is  a  function  of  the  Nyquist  bandlimit  of  the  segmented  wavefront 
control  device.  This  finding  leads  to  an  additional  goal  to  identify  the  relationship 
between  turbulence  strength  and  the  spatial  frequency  of  a  device. 
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1.3  Thesis  Outline 


This  thesis  meets  these  goals  through  a  detailed  investigation.  Chapter  II  ex¬ 
amines  the  basic  theoretical  properties  of  the  atmosphere  and  how  phase  screens  are 
used  to  represent  atmospheric  turbulence.  It  also  discusses  how  different  segmented 
wavefront  control  devices  alter  the  phase  of  incident  light  and  emulate  turbulence. 
Additionally,  Chapter  II  includes  a  detailed  literature  review  that  outlines  limitations 
others  have  found  for  segmented  wavefront  control  devices.  Chapter  III  then  looks  at 
the  laboratory  setup  of  wavefront  control  devices  to  determine  the  best  placement  of 
these  devices  to  maximize  the  range  of  turbulence  conditions  they  can  emulate.  Chap¬ 
ter  III  also  presents  a  theoretical  derivation  of  a  numerical  solution  for  the  maximum 
possible  amount  of  turbulence  as  a  function  of  the  number  of  pixels  on  a  segmented 
wavefront  control  device.  Chapter  IV  presents  a  complementary  analysis  using  com¬ 
puter  simulation  with  the  Matlab®  computer  program  from  The  MathWorks  Inc. 
Finally,  Chapter  V  gives  a  comparison  between  the  theoretical  and  simulation  results. 
This  chapter  then  provides  concluding  remarks  as  well  as  possible  extensions  of  this 
work  in  future  research. 
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II.  Theory  and  Review 


This  work  builds  on  that  of  previous  researchers,  moreover,  this  section  reviews 
what  others  have  done  in  the  field  of  atmospheric  propagation  using  SLMs.  It 
begins  with  basic  theory  concerning  atmospheric  turbulence  and  then  explores  other 
important  and  relevent  areas  that  are  the  foundations  of  this  research. 

2.1  Atmospheric  Properties 

Light  traveling  in  a  vacuum  can  easily  be  modeled  using  the  Fresnel  propagation 
integral.  [13]  However,  when  light  travels  through  the  Earth’s  atmosphere,  it  is  affected 
by  the  turbulence  present.  Understanding  how  that  turbulence  alters  the  path  of 
light  traveling  through  it  has  been  the  subject  of  many  publications.  One  of  the 
simplest  ways  to  understand  the  properties  of  atmospheric  turbulence  is  to  view  the 
refractive  index  of  the  atmosphere  as  a  random  process.  [1]  Random  variations  in  the 
temperature  and  pressure  of  the  Earth’s  atmosphere  spatially  and  temporally  alter 
the  atmosphere’s  index  of  refraction  n  for  a  given  point  R  in  space  according  to 

n{ R)  =  1  +  77.6  x  1(T6(1  +  7.52  x  10~3A“2)^®,  (2.1) 

T(R) 

where  A  is  the  optical  wavelength  in  micrometers,  P  is  the  pressure  in  millibars,  and 
T  is  the  temperature  in  Kelvins.  [1]  It  is  impossible  to  know  the  refractive  index  of  the 
atmosphere  at  all  points  in  space  at  every  moment  in  time.  Therefore,  a  statistical 
analysis  is  needed  to  best  represent  the  atmosphere.  [33] 

This  analysis  is  based  on  utilizing  turbulent  eddies  (i.e.,  pockets  of  air  with 
highly  correlated  indices  of  refraction).  [31]  Turbulence  can  be  defined  as  randomly 
distributed  eddies  of  varying  sizes  and  temperatures.  How  light  travels  through  the 
changes  in  index  of  refraction  due  to  turbulent  eddies  in  its  path  defines  how  it  is 
affected  by  atmospheric  turbulence.  The  outer  scale  L$  is  the  average  size  of  the 
largest  eddies.  Likewise,  the  inner  scale  Z0  is  the  average  size  of  the  smallest  eddies. 
As  inertial  forces  act  on  the  atmosphere,  the  larger  eddies  break  up  into  smaller  ones 
and  a  continuous  distribution  of  eddy  sizes  can  be  found  between  the  outer  and  inner 


5 


scales.  This  range  of  eddy  sizes  is  known  as  the  inertial  subrange.  When  eddies  break 
down  smaller  than  Iq,  they  are  considered  too  small  to  follow  the  statistics  of  turbulent 
flow  and  turbulence  gives  way  to  laminar  flow.  [1,33] 

The  simplest  method  of  understanding  the  random  fluctuations  in  index  of  re¬ 
fraction  is  to  use  its  structure  function  Dn(R).  This  function  can  be  used  to  describe 
the  mean-square  difference  in  index  of  refraction  between  two  points  separated  by  a 
distance  R  in  space.  The  structure  function  is  defined  by 

-Dra(Ri,  R2)  =  ([n(Ri)  —  n(R2)]")  ,  (2-2) 

where  Ri  and  R2  are  vectors  describing  points  in  space  and  (•)  is  an  ensemble  av¬ 
erage.  [1]  Assuming  that  refractive  index  fluctuations  maintain  stationary  increments 
further  simplifies  this  equation.  If  the  field  is  statistically  homogenous,  the  covariance 
function  Rn(R)  given  by 

£n(Ri,  R2)  =  (ni(Ri)ni(Ri  +  R)) ,  (2-3) 

becomes  independent  of  spatial  position  and  is  only  dependent  on  the  distance  be¬ 
tween  the  two  points  of  interest:  R  =  Ri  —  R2.  If  the  random  held  is  also  isotropic, 
the  covariance  function  becomes  dependent  only  on  the  scalar  distance  R  =  |Ri  — R2|2 
between  the  two  points.  The  structure  function  is  related  to  the  covariance  function 
by 

Dn{R)  =  2[Bn(0)  -  Bn(R)].  (2.4) 

When  these  assumptions  are  combined  with  Kolmogorov’s  dimensional  analysis,  the 
structure  function  is  reduced  to 

Dn(R)  =  C2nR2/z  for  l0  <  R  <  L0,  (2.5) 
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where  C2  is  the  refractive  index  structure  constant  and  has  units  of  m-2/3.  This  pa¬ 
rameter  defines  the  strength  of  the  fluctuations  in  the  refractive  index  and  is  generally 
referenced  as  the  strength  of  the  turbulence.  [1,31] 


The  power  spectral  density  (PSD)  of  the  random  changes  in  the  atmosphere’s 
index  of  refraction  <£„(#«)  also  characterizes  the  statistical  distribution  of  the  size  and 
number  of  turbulent  eddies  along  the  optical  path.  The  Wiener-Khintchine  theorem, 
given  by 

/OO 

Bn(R)e~j2™RdR,  (2.6) 

-OO 


states  that  the  power  spectrum  is  the  Fourier  transform  of  the  covariance  of  the  refrac¬ 
tive  index  variation  from  Eq.  (2.3).  [1]  This  theorem  and  Eq.  (2.4)  can  be  combined, 
leading  to  an  expression  for  the  refractive  index  structure  function  in  terms  of  the 
PSD,  given  by 


Dn(R)  =  87 r  /  ft2<Pn(ft) 


sin  (/cP) 

ft/? 


dtz, 


(2,7) 


where  ft  is  the  radial  polar  component  of  the  three-dimensional  angular  spatial  fre¬ 
quency  vector  ft),  measured  in  rad/m.  The  2/3  power-law  shown  in  Eq.  (2.5)  substi¬ 
tuted  into  Eq.  (2.7)  leads  to 


$„(k)  =  0.033C^ft-11/3  for  1/L0  <C  ft  <  l/l0.  (2.8) 


Eq.  (2.8)  is  more  commonly  known  as  the  Kolmogorov  power-law  spctrum  and  is  only 
valid  over  the  inertial  subrange.  Other  spectra  were  introduced  in  order  to  increase 
the  valid  range  for  the  PSD.  These  include  the  von  Karman,  Tatarskii,  and  modified 
spectrum  among  others.  The  modified  von  Karman  spectrum  is  given  by 

for  0  <  ft  <  oo, 

where  fto  =  27t/Lo  and  ftm  =  5.92/Zo-  [1]  Four  common  PSDs  are  shown  as  a  function 
of  spatial  frequency  in  Fig.  2.1.  For  simplicity  and  proof  of  concept,  this  research  will 


(2.9) 


4>„(k)  =  0.033C2  (ft  +  ft) 


-11/6 


exp 
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Figure  2.1:  Common  atmospheric  refractive  index  PSD  models  given  as  a  function  of 

spatial  frequency.  This  image  is  reproduced  with  the  consent  of  Schmidt  [33] 


use  the  Kolmogorov  power-law  spectrum  from  Eq.  (2.8)  to  represent  the  statistical 
properties  of  the  atmosphere. 

Other  important  atmospheric  parameters  are  the  Fried  parameter  r0  and  the 
Rytov  number  which  are  defined  by 


r0 


a 


2 

v. 


0.423&2  f  Cl(z)dz 
Jo 


-3/5 


0.5631  k7/6 


-  zf^dz, 


(2.10) 

(2.11) 


where  A  is  the  optical  wavelength  and  k  =  2tt/\.  Equations  (2.10)  and  (2.11)  are 
both  valid  for  a  plane  wave  propagating  from  z  =  0  to  z  =  L,  where  L  is  the  optical 
path  length  measured  in  meters.  [1]  The  Fried  parameter  defines  the  roll  off  of  the 
volume  underneath  the  optical  transfer  function  and  is  a  measure  of  the  coherence 
width  of  the  atmosphere.  The  Rytov  number  is  a  measure  of  the  scintillation  in  light 


as  it  travels  through  the  atmosphere.  [1,12]  It  is  equivalent  to  log-amplitude  variance 
under  weak  turbulence  conditions.  The  Fried  parameter  and  Rytov  number  are  direct 
measures  of  the  phase  and  amplitude  properties  of  an  optical  wave. 

2.2  Modeling  the  Atmosphere 

2.2.1  Layered  Model  of  Turbulence. 

Modeling  the  atmosphere  would  be  easy  if  it  were  possible  to  use  a  simple 
statistical  model.  Unfortuantely  this  is  not  usually  an  option  because  the  equations 
are  often  analytically  intractable  for  the  most  general  theories.  This  is  especially  true 
in  the  case  of  modeling  complex  AO  systems.  [31]  The  most  common  solution  to  this 
problem  is  to  break  the  turbulent  optical  path  into  a  finite  number  of  discrete  layers. 
Under  certain  constraints,  a  layered  model  is  valid  for  an  incident  plane  wave.  [31] 

Each  layer  is  a  thin  phase  screen  that  represents  a  much  thicker  volume  of  at¬ 
mospheric  turbulence.  As  light  travels  through  a  phase  screen,  it  accumulates  the 
turbulence  statistics  associated  with  a  long  propagation  through  turbulence.  The 
phase  screen  must  be  much  thinner  than  the  section  of  the  total  propagation  distance 
it  represents.  To  propagate  a  wavefront  through  a  simulated  turbulence  path  using 
phase  screens,  the  light  is  first  propagated  from  the  source  to  the  first  screen  as  if 
it  were  traveling  through  a  vacuum.  Next,  the  complex  transmittance  function  of 
the  first  phase  screen  is  multiplied  to  the  incident  wavefront.  The  wavefront  contin¬ 
ues  propagating  to  the  next  phase  screen  as  if  through  a  vacuum,  and  this  process 
repeats  until  the  wavefront  reaches  the  target.  This  is  known  as  the  split-step  beam- 
propagation  method  and  is  shown  for  a  two  phase  screen  system  in  Fig.  2.2.  [33] 

The  split-step  beam-propagation  approach  allows  for  simulation  of  a  broad  and 
continuous  atmosphere  in  a  discrete  method.  The  ith  layer  in  the  model  can  be  the¬ 
oretically  defined  by  the  structure  parameter  C%.,  the  height  above  the  ground  Zj, 
and  the  thickness  A z,  of  the  layer  of  atmospheric  turbulence  it  represents.  Values 
for  these  parameters  are  chosen  so  that  the  low-order  statistical  moments  m  for  the 


9 


Phase  Phase 


Figure  2.2:  Split-step  beam  propagation  method  using  two  phase  screens.  Plane-wave 

light  is  incident  from  the  left  and  propagates  through  vacuum  to  the  first  phase  screen. 
The  light  is  distorted  by  the  first  screen  and  then  propagates  through  vacuum  to  the  second 
screen.  The  second  screen  further  distorts  the  light  and  it  then  travels  through  vacuum  to  the 
observation  plane.  The  resultant  wavefront  at  the  observation  plane  has  the  same  turbulence 
characteristics  as  light  that  has  traveled  continuously  through  a  turbulent  atmosphere. 


continuous  model  match  those  for  the  layered  model 

rL  N 

/  C2n(z')(zTdz'  =  '£Clz?  A*,  (2.12) 

i= 1 

where  L  is  the  total  propagation  distance  and  N  is  the  number  of  phase  screens 
used  to  represent  L.  [31]  This  allows  for  important  atmospheric  parameters  to  be 
computed  with  a  discrete  sum.  [33]  The  Fried  parameter  r0  and  Rytov  number  cr\  are 
represented  continuously  by  Eqs.  (2.10)  and  (2.11),  respectively,  and  can  be  computed 
for  a  layered  turbulence  model  using 


r0 


a 


2 

n 


0A23k2J2c2niAZi 

i 


-3/5 


0.5631  k7/6J2CUL  ~  -*)5/6A^- 

i 


(2.13) 

(2.14) 
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An  effective  coherence  diameter  r0i  can  be  found  for  the  ith  turbulent  layer  by  com¬ 
bining  terms  in  Eq.  (2.13)  and  is  given  by 

r0i  =  [0A23k2Cl.Azi]  "3/5  .  (2.15) 

The  coherence  diameter  for  each  layer  can  then  be  combined  to  find  the  overall  Fried 
parameter  for  the  entire  length  of  atmosphere  being  represented  and  is  given  by 

/  N  \  -3/5 

r°=  (Er»7/3j  •  <2-16) 

The  Fried  parameter  is  often  used  in  addition  to  the  structure  function  to  define  the 
strength  of  the  turbulence  in  a  phase  screen.  [33]  Using  a  similar  analysis,  the  overall 
Rytov  number,  given  by 

N 

4  =  1-331V5/6  ]T  r-5/s(i  -  (2.17) 

2=1 

can  also  be  written  in  terms  of  the  individual  coherence  diameters  for  each  layer.  As 
also  seen  in  the  continuous  model  defined  in  Section  2.1,  it  should  be  noted  that  the 
above  equations  are  only  valid  for  a  propagating  plane  wave  traveling  from  z  =  0  to 
2  =  L.  [1] 

2.2.2  Phase  Screen  Generation:  Fourier  Series  Method. 

Phase  screens  are  created  using  a  computer  generated  two-dimensional  array  of 
random  numbers  that  have  been  manipulated  into  a  grid  of  phase  values  with  the 
same  properties  as  the  atmosphere  they  represent.  The  method  used  in  this  research 
to  manipulate  those  random  values  and  create  phase  screens  is  the  Fourier  series 
approach. 

The  first  step  to  creating  phase  screens  using  the  Fourier  series  method  is  to 
choose  the  PSD  for  the  turbulence  desired  based  on  L0,  /0,  r0,  and  the  linear  spatial 
frequency  (fx,  fy).  [7]  The  PSD  must  be  sampled  at  a  fast  rate  in  the  low  frequency 
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region  and  less  often  in  the  high  frequency  region  in  order  to  adequately  sample  the 
spectrum.  This  gives  more  weight  to  the  lower  order  aberrations  that  are  responsible 
for  most  of  the  turbulence.  [21] 

Next,  Fourier  series  coefficients  are  found  for  the  spatial  frequency  region  of 
interest.  This  is  done  by  randomly  drawing  Gaussian  variables  with  the  appropriate 
variance  given  by  the  PSD.  The  Fourier  series  coefficients  are  denoted  by  c„)TO.  Each 
has  circular  Gaussian  statistics  with  variance  corresponding  to  the  sampled  PSD. 
Each  randomly  drawn  array  of  the  complex  coefficients  represents  a  unique  random 
phase.  A  Fourier  series  phase  (pk{x,y)  is  calculated  by 

OO  OO 

Mx,y)=  E  cn,mej2n{f™x+f^y\  (2.18) 

n= — oo  m=—oo 

where  (x,y)  are  the  spatial  coordinates  of  the  phase  screen  and  ( fXn,fym )  are  the 
spatial  frequency  components.  One  of  the  main  benefits  to  this  approach  is  that  the 
phase  screen  is  defined  for  all  points  in  space.  The  coefficients  have  to  be  deter¬ 
mined  only  once  and  then  using  the  spatial  frequency  limits,  the  phase  screen  can  be 
synthesized  anywhere  in  space.  [21] 

2.3  Zernike  Polynomials 

While  any  basis  (e.g.  Fourier  series)  can  be  used  to  represent  aberrations  in 
a  phase  screen,  the  Zernike  series  is  selected  for  this  research.  This  is  because  it 
is  orthogonal  over  a  circular  aperture.  Random  aberrations  that  result  from  light 
traveling  through  turbulence  can  be  represented  by  a  weighted  sum  of  Zernike  poly¬ 
nomials.  [31]  Zernike  modes  are  two-dimensional  polynomials  that  are  represented 
here  as  orthonormal  on  an  aperture  of  diameter  D.  The  Zernike  polynomials  Zmjn 
are  given  in  polar  coordinates  by 

Zm.,n(P->e)x  =  Vn  +  l  R™(2p/D)V2cos(m6)  (2.19) 

Zm,n(p,6)y  =  VnTl  Rn{2p/D)V2sin(m9)  (2.20) 
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Table  2.1:  First  eleven  Zernike  polynomials  and  their  common  names. 


Mode  (z) 

m 

n 

Zm,n(pi 

Name 

1 

0 

0 

l 

piston 

2 

1 

1 

2(2 p/D)  cos(P) 

x  tilt 

3 

1 

1 

2(2 p/D)  sin(6)) 

y  tilt 

4 

0 

2 

y/3  [2(2 p/D)2  -  1] 

defocus 

5 

2 

2 

v/6(2p/D)2cos(20) 

x  primary  astigmatism 

6 

2 

2 

VQ(2p/D)2sm{29) 

y  primary  astigmatism 

7 

1 

3 

V8  [3(2 p/D)3  -  2(2 p/D)}  cos (9) 

x  primary  coma 

8 

1 

3 

V8  [3(2 p/D)3  -  2(2 p/D)}  sin(0) 

y  primary  coma 

9 

3 

3 

y/8(2p/D )3  cos  (30) 

x  trefoil 

10 

3 

3 

V/8(2p/H)3sin(36*) 

y  trefoil 

11 

0 

4 

y/5  [6(2p/ D)4  -  6(2 p/D)2  +  1] 

primary  spherical 

for  m^O,  and 

Zo,n(P,  e)  =  R°n(2p/D)  (2.21) 

for  m  =  0,  where  p  is  the  transverse  component  of  the  two-dimensional  spatial  vector 
R.  The  radial  functions  R™(2p/D)  used  in  these  equations  are  defined  by 


K(?p/d)  = 


(n—m)/ 2 

E 

q= o 


-l)5 


(n  —  q)\(2p/D)n  2q 


q\  [[n  +  m)/2  —  q\\ [[n  —  m)/2  —  q\\ ’ 


(2.22) 


where  the  azimuthal  order  m  and  radial  order  n  are  non-negative  integers  that  satisfy 
m  <  n  and  n  —  m  =  even.  [31,32]  These  polynomials  can  be  found  for  any  combination 
of  (m,n),  however,  the  low-order  modes  represent  the  most  common  aberrations  that 
affect  optical  systems.  Zernike  modes  are  sometimes  referenced  by  their  mode  index 
number  i.  Some  common  low-order  Zernike  modes  are  given  in  Table  2.1,  and  a 
visual  example  of  some  common  modes  present  in  a  turbulent  phase  screen  is  given 
in  Fig.  2.3. 
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(d)  (e)  (f) 

Figure  2.3:  Examples  of  Zernike  aberrations  present  in  a  turbulent  wavefront  with  the 

following  modes  present:  (a)  all  Zernike  modes  (b)  piston  (c)  X  and  Y  tilt  (d)  X  and  Y 
astigmatism  (e)  X  and  Y  coma  (f)  tertiary  spherical. 


The  phase  of  a  wavefront  4>(Dp,  6)  can  be  represented  as  a  sum  of  all  the  Zernike 
modes  in  the  wave  with 


N 


<t>{Dp,Q)  =  ^aiZi(p,6),  (2.23) 

2=1 

where  a*  is  the  weight  associated  with  the  ith  Zernike  polynomial  and  N  is  the  number 
of  summed  Zernike  modes.  The  expression  used  to  find  the  weights  for  individual 
Zernike  modes  in  cartesian  coordinates  within  a  given  wavefront  y )  is 


ffZ  W (x’  y)zi(x ,  y)dxdy 

f/Z  W (x’  y)ZKx,  y)dxdy 


(2.24) 


where  W (x,  y)  is  the  aperture  function  for  a  circle  of  unit  radius.  [31] 

Equation  (2.23)  is  verified  in  Fig.  2.4  which  shows  that  as  N  increases,  the  phase 
of  the  wavefront  approaches  the  complete  turbulent  phase  with  all  modes  included. 
Figure  2.4(f)  is  the  turbulent  wavefront  created  by  summing  the  first  37  Zernike 
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(d)  (e)  (f) 

Figure  2.4:  Turbulent  wavefront  phase  with  the  following  Zernike  modes  present:  (a)  all 
Zernike  modes  (b)  piston  only  (c)  piston  through  X  and  Y  tilt  (d)  piston  through  X  and 
Y  astigmatism  (e)  piston  through  X  and  Y  coma  (f)  piston  through  tertiary  spherical.  All 
subfigures  are  plotted  on  the  same  scale. 

modes.  It  can  be  seen  that  Fig.  2.4(f)  begins  to  resemble  Fig.  2.4(a)  in  which  all 
Zernike  modes  are  present.  Similarly,  Eq.  (2.23)  can  be  used  to  show  that  subtract¬ 
ing  Zernike  modes  from  a  wavefront  leaves  only  the  residual  phase.  As  that  phase 
approaches  to  zero,  the  wavefront  loses  its  turbulence  and  becomes  a  plane  wave. 
Figure  2.5  presents  an  example  of  this  phenomenon.  It  can  be  seen  that  Fig.  2.5(f) 
is  approaching  a  non-aberrated  wavefront  because  it  is  has  lost  nearly  all  of  its  phase 
and  it  resembles  an  AO  compensated  wavefront. 

To  synthesize  a  phase  screen  using  the  Zernike  series,  the  coefficients  are  com¬ 
puted  beginning  with  random  draws  from  the  standard  normal  probability  density 
function.  These  are  then  transformed  to  have  the  covariance  given  by  Eqs.  (2.25)  and  (2.26). 
An  expression  derived  from  the  Kolmogorov  power  spectrum  for  the  covariance  be- 
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Figure  2.5:  Residual  phase  of  a  turbulent  wavefront  with  the  following  Zernike  modes 

removed:  (a)  none  (b)  piston  only  (c)  piston  through  X  and  Y  tilt  (d)  piston  through  X 
and  Y  astigmatism  (e)  piston  through  X  and  Y  coma  (f)  piston  through  tertiary  spherical. 
All  subfigures  are  plotted  on  the  same  scale. 


tween  the  Zernike  coefficients  is  given  by 


(oiOj)  =  0.0072 


5/3 


X 


(_1)K+nj_2mi)/2  [(n.  +  1)(n  .  +  1}]l/2  ^8/3^ 

_ r(14/3)r  [{rij  +  rij  —  5/3)/2] _ 

r  [(m  -  rij  + 17/3)/2]  r  [(nj  -  m  + 17/3)/2]  r  [(m  +  nj  +  23/3) /2] 


,  (2.25) 


for  i  —  j  =  even,  and 


(0*%)  =  0 


(2.26) 


for  i—  j  —  odd.  [25,31]  In  Eq.  (2.25)  Smimj  is  the  Kronecker  delta  function  and  T(x) 
is  the  Gamma  function  given  by  [1] 


r(aO 


e-Hx-1 


dt , 


Re(x)  >  0. 


(2.27) 
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An  approximation  for  the  mean  square  residual  phase  error  e2  when  modes  1  through 
i  are  removed,  was  found  by  Noll  [25]  and  is  given  in  rad2  by 

e2  «  0.29Mi-^/2(D/ro)5/3.  (2.28) 

Equation  (2.28)  is  only  valid  for  higher  order  modes  {i  >  10),  so  in  his  article,  Noll  [25] 
presents  the  results  from  a  laborious  technique  performed  by  Fried  to  calculate  the 
residual  phase  found  in  the  lower  order  modes. 

2.4  Wavefront  Control  Devices 

Thus  far,  wavefront  control  has  only  been  discussed  abstractly.  In  practice, 
there  are  two  ways  that  devices  alter  the  phase  of  a  wavefront.  They  can  either 
alter  the  index  of  refraction  that  the  wavefront  travels  through  or  change  the  distance 
traveled  by  the  light.  LC  SLMs  use  the  property  of  birefringence  to  alter  the  index 
of  refraction  in  the  SLM.  [13]  In  contrast,  DMs  use  tiny  actuators  to  alter  the  shape 
of  the  wavefront  so  that  the  physical  distance  traveled  by  light  is  modified. 

2-4-1  Liquid- crystal  Spatial  Light  Modulators. 

Atmospheric  turbulence  is  characterized  by  the  change  in  index  of  refraction 
of  the  atmosphere.  LC  materials  have  been  utilized  in  order  to  duplicate  this  effect 
and  create  phase  screens.  LC  devices  are  commonly  used  in  laboratory  environments 
for  adaptive  optics  because  they  are  nonmechanical,  inexpensive,  and  highly  reliable. 
LC  materials  share  properties  of  both  liquids  and  solids.  These  elliptically-shaped 
molecules  that  compose  LC  materials  are  not  bound  to  each  other  in  a  rigid  formation, 
so  they  are  able  to  rotate  or  slide  when  an  electrical  or  mechanical  force  is  applied. 
LC  materials  come  in  a  variety  of  meso-phases;  however,  the  one  of  interest  to  this 
research  is  the  nematic  liquid  crystal  (NLC).  The  molecules  in  a  NLC  tend  to  align 
in  parallel  with  their  centers  randomly  distributed  throughout  the  crystal.  [13] 
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The  molecules  in  these  devices  have  positive  uniaxial  birefringent  characteristics 
that  rotate  when  an  electric  field  is  applied.  Light  traveling  through  the  elliptical 
molecules  of  a  NLC  SLM  experiences  an  altered  index  of  refraction  based  on  the 
alignment  of  the  molecules.  The  NLC  SLM  is  composed  of  several  layers.  Linearly 
polarized  incident  light  first  passes  through  a  transparent  glass  cover,  then  through  a 
transparent  electrode,  and  then  into  the  NLC  material.  The  incident  light  obtains  a 
phase  delay  as  it  travels  through  the  NLC  material  and  then  reflects  from  an  array  of 
electrodes  before  traveling  back  out  along  the  same  path.  The  double  pass  through 
the  NLC  material  causes  a  phase  delay  corresponding  to  the  electric  field  applied  to 
that  area  of  the  device  and  the  reflective  electrode  array  subdivides  the  NLC  SLM  into 
an  array  of  independently  controlled  elements.  Each  element  can  impart  a  different 
phase  delay  on  the  incident  wavefront.  Large  XY  arrays  of  control  circuitry  behind 
the  NLC  SLM  allow  for  hundreds  of  thousands  of  controllable  phase  alterations  in 
current  commercial  devices.  [13,33]  Figure  2.6  provides  an  example  of  how  an  applied 
voltage  alters  the  molecular  alignment  inside  a  transmissive  NLC  SLM.  [33] 

2-4-1 -1  Practical  issues  for  adaptive  optics. 

The  ability  of  a  NLC  SLM  to  impart  many  controllable  phase  distortions  on  a 
single  wavefront  makes  it  useable  as  a  phase  screen.  In  1998,  Dayton  et  al.  [9]  were 
among  the  first  to  use  a  NLC  SLM  as  a  turbulent  phase  screen.  That  group  used  a 
NLC  device  created  by  Meadowlark  Optics  Inc.  with  an  11  mm  active  area  composed 
of  127  cells  arranged  hexagonally.  It  was  found  that  the  device  was  capable  of  produc¬ 
ing  up  to  36  Zernike  modes  and  was  thus  capable  of  providing  good  approximations 
for  Kolmogorov  phase  fronts  with  D/r0  <  10.  These  measurements  were  analyzed 
using  the  structure  function  of  the  altered  wavefront  and  the  error  in  Zernike  mode 
representation.  It  was  found  that  spatial  quantization  limited  the  order  of  spatial  cor¬ 
rection  that  could  be  achieved.  In  addition  to  Zernike  mode  representation,  Dayton 
et  al.  showed  the  ability  of  a  NLC  SLM  to  compensate  time  evolving  phase  screens 
as  an  adaptive  corrector.  One  year  earlier,  Love  [20]  was  able  to  produce  15  Zernike 
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Figure  2.6:  Diagram  of  using  an  applied  voltage  to  control  the  birefringence  inside  a  NLC 
SLM.  The  molecules  rotate  with  respect  to  the  applied  voltage  on  the  device.  Incident  light 
travels  from  the  top  and  passes  through  the  transparent  electrode  before  hitting  the  rotated 
NLC  molecules.  The  rotated  molecules  in  the  middle  of  the  schematic  cause  a  phase  delay 
in  the  transmitted  light  with  respect  to  the  edges  because  the  material  is  uniaxial  positive. 

This  image  is  reproduced  with  the  consent  of  Schmidt  [33] 


modes  on  a  69  pixel  device  that  was  also  made  by  Meadowlark  Optics  Inc.  His  work 
proved  that  quality  optical  devices  could  be  made  using  NLC  SLMs.  It  is  common 
for  today’s  NLC  SLMs  to  have  at  least  512  x  512  pixels. 

Brooks  [7]  and  Phillips  [27]  further  demonstrated  the  abilities  of  NLC  SLMs. 
Both  researchers  used  a  512  x  512  pixel  NLC  SLM  created  by  Boulder  Nonlinear 
Systems  Inc.  with  a  7.68  mm  active  area.  Brooks  was  able  to  demonstrate  that 
atmosphere-like  aberrations  could  be  added  to  an  optical  system  via  a  NLC  SLM. 
He  was  one  of  the  first  to  perform  atmospheric  simulation  using  one  of  these  devices. 
Phillips  was  able  to  determine  that  a  system  could  be  created  using  multiple  SLMs  to 
independently  control  the  Fried  parameter,  isoplanatic  angle,  and  the  Rytov  variance 
of  the  simulated  atmosphere.  Both  researchers  demonstrated  that  experimental  values 
could  be  scaled  appropriately  to  show  good  correlation  with  the  theoretical  and 
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simulation  values.  Brooks  created  a  probability  density  function  of  irradiance  from 
all  phase  screen  realizations  and  compared  it  to  the  theoretically  expected  log-normal 
distribution  to  verify  the  match.  He  found  that  even  though  the  error  between  the 
experimental  and  theoretical  cases  was  fairly  high,  the  intensity  distributions  still  had 
the  form  and  properties  of  log-normal  curves.  Phillips  also  performed  an  intensity 
evaluation  to  verify  that  the  Rytov  variance  could  be  adequately  modeled  using  only 
two  phase  screens.  He  found  that  the  system  creates  higher  Rytov  variances  than 
expected  for  weak  turbulence  and  lower  variations  than  expected  for  the  opposite 
case.  The  error  in  the  experimental  simulation  ranged  from  2%  for  nominal  turbulence 
to  75%  for  very  weak  turbulence.  This  was  attributed  to  the  SLM  being  driven  to 
its  limits  (in  terms  of  applied  voltage)  for  the  very  strong  and  very  weak  turbulence 
cases,  thereby  causing  significant  error.  Brooks  also  performed  a  phase  evaluation 
to  further  validate  the  abilities  of  the  NLC  SLM.  The  phase  of  the  wavefront  in  the 
observation  plane  was  compared  to  the  average  structure  function  of  each  scenario 
over  many  realizations  with  the  analytic  Kolmogorov  structure  function.  Brooks  found 
experimental  phase  measurement  errors  from  15.7-26.6%  which  were  determined  to 
be  average  quality  results.  Despite  the  large  error,  it  was  shown  that  the  experimental 
structure  function  still  adhered  to  the  5/3  power  law.  [7,27] 

One  of  the  most  practical  issues  for  NLC  SLMs  in  adaptive  optics  is  the  prob¬ 
lems  that  can  come  from  not  properly  calibrating  the  device.  Both  Brooks  [7]  and 
Phillips  [27]  mentioned  that  diffraction  from  the  square  apertures  of  the  SLMs  may 
have  been  a  significant  contributing  factor  to  the  errors  seen  in  the  intensity  and  phase 
evaluations.  However,  poor  calibration  of  the  LC  SLMs  was  not  presented  as  a  large 
source  of  error  in  either  paper.  In  his  paper,  Schmidt  [35]  discusses  how  improper 
calibrations  can  lead  to  significant  error. 

Schmidt  used  the  same  LC  SLM  as  Brooks  [7]  and  Phillips  [27].  This  time,  it 
was  used  to  show  the  effects  of  proper  calibration  on  the  wavefront  in  the  observation 
plane.  He  showed  that  a  phase-to-command  calibration  method  could  be  improved 
upon  by  introducing  a  quarter-wave  plate  into  the  optical  path  and  performing  a 
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two-step  calibration.  First,  the  irradiance  was  measured  without  the  wave  plate,  and 
then  it  was  measured  with  the  wave  plate.  This  extra  step  preserved  the  sign  of 
the  derivative  of  the  measured  phase  and  allowed  for  a  standard  phase  unwrapping 
program  to  compute  the  true  physical  phase  written  onto  the  SLM.  In  addition  to 
improving  the  method  of  calibration  of  the  phase  loaded  onto  each  pixel  of  the  SLM, 
this  was  the  first  time  a  theoretical  analysis  of  the  amplitude  modulation  phase-to- 
command-value  calibration  was  shown.  [35] 

Schmidt  [35]  also  developed  a  method  of  calibrating  the  SLM  for  the  static 
aberration  inherent  in  every  device  due  to  a  slightly  warped  backplane  formed  dur¬ 
ing  the  manufacturing  process.  Most  SLMs  are  found  to  have  one  to  three  peak 
to  valley  waves  of  static  aberration.  The  static  aberration  in  the  SLM  described  in 
Schmidt’s  [35]  paper  was  measured  using  a  Twyman-Green  interferometer  and  com¬ 
pensated  by  commanding  a  nominal  phase  map  onto  the  device.  The  true  benefits 
of  having  a  properly  calibrated  device  were  then  demonstrated  in  the  paper.  [35]  In 
another  article,  Hart  [16]  even  showed  that  static  aberrations  on  SLMs  can  be  char¬ 
acterized  and  removed  using  simple  and  widely  available  hardware  without  requiring 
expensive  aberration  sensing  equipment.  Any  experiment  using  a  SLM  should  ensure 
proper  calibration  before  analyzing  any  results. 

2-4-2  Noted  Limitations. 

Several  limitations  exist  in  the  use  of  NLC  SLMs  as  adaptive  optics.  The  NLC 
material  itself  is  wavelength  dependent.  Therefore,  different  incoming  wavefronts  will 
see  different  indices  of  refraction  and  optical  path  lengths  when  they  hit  the  device. 
NLC  SLMs  are  manufactured  to  operate  within  specific  wavelength  regions.  These 
devices  also  only  operate  correctly  with  linearly  polarized  light.  Bold  et  al.  [6]  address 
this  issue  with  a  discussion  of  how  a  SLM  can  be  coupled  with  a  1/4- wave  plate  to 
provide  uniform  modulation  of  unpolarized  light.  SLMs  also  suffer  a  decrease  in  signal 
throughput  and  an  increase  in  diffraction  artifacts  as  a  result  of  the  opaque  transistors 
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surrounding  each  pixel.  [3]  However,  these  effects  can  be  compensated  to  some  degree 
through  proper  calibration  of  the  device.  [16] 

The  most  prominent  issue  that  plagues  NLC  SLMs  is  the  relaxation  time  of  the 
material.  When  a  phase  is  electronically  commanded  to  the  device,  it  can  respond 
extremely  fast.  However,  the  molecules  in  the  NLC  material  take  comparatively 
much  longer  to  relax  back  to  their  original  state  after  that  command  is  removed. 
Current  devices  have  response  times  on  the  order  of  several  milliseconds.  Serati  et 
al.  [36]  present  theory  proving  the  possibility  of  sub-millisecond  response  times  for 
NLC  SLMs.  Although  close,  manufacturers  still  have  not  been  able  to  create  devices 
that  are  capable  of  such  frame  rates.  There  is  much  work  in  high-speed  LC  materials 
in  which  dual-frequency  devices  have  been  moderately  successful.  [8,24,29,36]  This 
limitation  constricts  SLMs  to  a  laboratory  environment  and  prevents  them  from  being 
used  in  real-time  high  speed  AO  applications. 

Litvin  et  al.  [19]  and  Bagwell  et  al.  [3]  also  discuss  how  digitization  is  another 
limiting  constraint  of  SLM  devices.  Most  current  SLMs  have  eight  bits  of  phase  per 
27t  of  phase.  The  desired  analog  phase  value  for  a  pixel  must  be  converted  to  a  number 
from  0  to  255  in  order  to  be  electronically  commanded  to  the  device.  The  polarity 
independence  of  NLC  molecular  rotation  further  limits  this  range  to  seven  bits  or 
128  programmable  levels.  This  restriction  exceeds  the  spatial  resolution  of  the  pixels 
on  most  devices  and  leads  to  sampling  as  the  limiting  constraint  when  using  NLC 
SLMs.  [3,19] 

2. 4- 2.1  Other  uses. 

Liquid-crystal  SLMs  can  be  used  in  a  vast  array  of  other  applications.  They  can 
be  used  for  AO  to  both  simulate  atmospheric  turbulence  and  correct  for  it.  It  has  also 
been  shown  that  SLMs  can  be  used  to  steer  a  laser  beam.  Recently,  Linncnberger  et 
al.  [18]  showed  that  LC  devices  can  be  used  to  steer  a  laser  beam  up  to  ±6.95  degrees. 
Hallstig  et  al.  [14, 15]  not  only  discusses  how  using  LC  SLMs  to  alter  the  phase  of  a 
wavefront  can  induce  a  deflection  that  acts  to  steer  a  laser  beam,  but  goes  on  to  show 
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that  a  SLM  can  be  used  to  alter  the  shape  of  the  beam  itself.  Schmidt  [34]  expands  on 
beam  shaping  through  the  use  of  a  phase  retrieval-based  wavefront  control  to  shape  a 
laser  beam  into  a  turbulence-distorted  beam.  This  effectively  eliminates  the  need  for 
multiple  phase  screens  in  a  layered  turbulence  model  and  allows  LC  SLMs  to  model 
bulk  atmospheric  turbulence  without  moving  parts.  Bagwell  et  al.  [3]  use  SLMs  to 
act  as  lenses  and  even  demonstrate  nonmechanical  zoom  and  enhanced  multi-spectral 
imaging  systems  using  liquid-crystal  devices. 

2-4-3  Deformable  Mirrors. 

The  DM  is  another  widely  used  device  to  correct  for  aberrations  in  a  wavefront. 
As  stated  previously,  DMs  use  actuators  located  behind  the  mirror  surface  to  alter 
the  shape  of  the  wavefront  so  that  the  physical  distance  traveled  by  light  is  modified. 
Changing  the  length  of  the  optical  path  creates  a  phase  delay  at  the  observation  plane. 
The  DM  acts  as  a  phase  screen  capable  of  commanding  many  phase  distortions  on 
a  single  wavefront  when  the  light  is  spread  over  many  actuators.  DM  actuators  are 
capable  of  providing  several  micrometers  of  stroke  with  nanometer  precision.  [4]  DMs 
are  the  preferred  device  used  in  real-time  astronomical  applications  due  to  their  fast 
operating  speed.  [11]  However,  such  DMs  are  often  too  expensive  to  use  in  a  laboratory 
environment.  Instead,  a  smaller  and  more  cost  efficient  device  was  needed  for  use  in 
research  applications.  MEMS  DMs  provide  a  solution  by  combining  technology  from 
the  fields  of  micro-instruments,  adaptive  optics,  and  controls  to  form  a  unique  mirror 
assembly.  [23]  These  mirrors  have  the  advantages  of  being  simple  to  manufacture, 
inexpensive,  lightweight,  and  integrable  with  drive  and  sensing  electronics.  [23]  They 
also  consume  very  little  power,  are  capable  of  being  driven  at  KHz  frame  rates,  and, 
unlike  SLMs,  most  DMs  are  wavelength-independent  due  to  the  reflective  properties 
of  the  mirror  surface.  [5, 11]  The  three  main  types  of  MEMS  DMs  in  use  today  are 
continuous  mirror,  segmented  mirror  with  piston  motion  only,  and  segmented  mirror 
with  tip/tilt  motion  DMs.  [5]  MEMS  DMs  commonly  have  5x5  to  12x12  pixels/mirror 


23 


segments  and  the  more  expensive  models  have  as  many  as  32  x  32  pixels/mirror 
segments. 


2. 4  -3.1  Continuous  DMs. 

Continuous  facesheet  membrane  DMs  are  created  by  stretching  a  thin  mirror 
membrane  over  a  circular  silicon  frame.  The  shape  of  the  membrane  is  then  altered 
through  the  use  of  electrostatic  forces.  The  first  continuous  membrane  MEMS  DM 
had  nine  actuators  and  was  designed,  manufactured,  and  tested  in  the  end  of  the  20th 
century  at  Boston  University.  [4].  Another  method  of  creating  a  continuous  facesheet 
DM  is  to  use  a  thin  glass  faceplate  that  is  controlled  with  piezoelectric  actuators. 
Xinetics  Inc.  is  one  company  that  manufactures  this  type  of  device. 

A  continuous  surface  offers  several  advantages  over  other  mirror  types.  The 
surfaces  of  both  the  membrane  and  thin  glass  continuous  facesheet  DMs  ensure  smooth 
and  continuous  phase  variations  over  the  entire  mirror.  A  smooth  surface  causes 
almost  no  diffraction  in  the  reflected  beam.  Also,  there  is  no  loss  of  optical  intensity 
due  to  fill  factors  present  in  an  array  like  the  NLC  SLM.  These  advantages  do  not 
come  without  drawbacks  to  the  continuous  DM.  For  example,  the  continuous  surface 
causes  interactuator  coupling  that  complicates  adaptive  optic  control  between  the 
actuators.  This  coupling  is  known  as  the  influence  function  and  is  a  compromise  of 
the  system.  An  influence  function  of  approximately  ten  percent  is  typical  between 
adjacent  actuators.  [5]  This  influence  function  causes  a  tension  or  compression  on  the 
mirror  surface  that  can  cause  adverse  affects  in  MEMS  devices. 

Fernandez  and  Artal  [11]  were  able  to  use  a  37  element  continuous  MEMS 
DM  manufactured  by  OKO  Technologies  as  an  adaptive  optical  corrector.  They 
successfully  reproduced  the  first  21  Zernike  modes  with  good  accuracy  using  the 
mirror.  Unfortunately,  the  limited  range  of  control  voltages  allowed  some  of  the  higher 
order  modes  to  only  be  producible  over  a  very  small  range.  Despite  this  limitation, 
it  was  a  successful  demonstration  of  adaptive  wavefront  control  using  a  continuous 
facesheet  MEMS  DM.  [11]  Doble  and  Williams  [10]  are  working  to  apply  continuous 
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MEMS  DMs  in  the  field  of  optometry.  They  discuss  how  glasses  and  contact  lenses 
are  static  devices  that  can  only  correct  low-order  aberrations  such  as  defocus  and 
astigmatism.  Doble  and  Williams  present  the  idea  that  continuous  MEMS  devices 
can  be  used  as  adaptive  optics  for  real-time  complete  correction  of  human  vision. 
In  order  to  make  this  concept  a  reality,  it  is  necessary  to  continue  to  improve  the 
MEMS  devices  that  are  commercially  available.  The  challenge  lies  in  manufacturing 
a  compact  and  reliable  DM  with  large  stroke  and  actuator  count  at  a  low  cost.  [10] 

2. 4 .3. 2  Segmented  DMs. 

Another  type  of  MEMS  DM  is  the  segmented  DM.  This  type  of  mirror  comes  in 
two  different  versions.  The  first  type  is  the  segmented  MEMS  DM  with  piston-only 
motion.  This  DM  is  created  with  tiny  mirror  segments  at  the  end  of  each  actuator 
which  are  completely  independent  of  each  other.  This  mirror  has  no  interactuator 
coupling  between  the  mirrors,  however  undesirable  diffraction  effects  are  introduced 
due  to  the  gaps  between  adjacent  mirror  segments.  Also,  this  mirror  does  not  allow 
for  spatially  continuous  phase  modulation  across  the  surface  of  the  mirror.  [5]  The 
benefit  to  this  mirror  is  that  it  allows  for  relief  of  the  interactuator  stress  prominent 
in  continuous  facesheet  DMs.  Furthermore,  it  also  provides  the  most  mirror  segments 
over  any  device  due  to  the  one-to-one  ratio  of  segments  to  actuators. 

The  second  type  of  segmented  MEMS  DM  is  a  hybrid  between  the  previous 
alternatives.  A  segmented  DM  with  tip/tilt  motion  has  mirror  segments  over  three 
actuators.  Sometimes  these  actuators  are  shared  between  segments.  In  either  case, 
this  device  allows  for  the  optical  phase  to  be  closely  matched  at  the  interface  between 
adjacent  segments.  [5]  There  are  some  diffraction  and  interactuator  coupling  issues 
associated  with  this  device,  but  each  is  an  acceptable  compromise  between  the  other 
two  types  of  MEMS  DMs.  Perreault  et  al.  [26]  were  able  to  use  segmented  DMs 
with  tip/tilt  motion  to  achieve  phase  modulation  in  an  AO  system.  Their  research 
showed  significant  reduction  in  the  wavefront  phase  error  after  correction.  The  Strehl 
ratio  of  the  wavefront  increased  from  0.0034  to  0.1950  when  the  aberrated  wavefront 
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was  compensated.  [26]  This  Strehl  ratio  is  a  performance  metric  used  in  imaging. 
Strelil  ratios  close  to  one  represent  a  near-diffraction  limited  beam  with  almost  no 
adverse  turbulence  affects.  The  segmented  mirrors  relieved  problems  caused  by  stress 
on  a  continuous  membrane  and  still  maintained  a  fill  factor  of  98.6%.  This  type  of 
mirror  also  maintained  phase  continuity  between  segments  through  actuator  sharing 
of  adjacent  mirror  segments.  [26] 

Lee  et  al.  [17]  were  able  to  show  correction  near  the  diffraction  limit  for  spher¬ 
ical  aberration  using  a  segmented  MEMS  DM.  The  group  corrected  an  aberrated 
wavefront  with  a  0.7  meter  radius  of  curvature  within  two  percent  of  the  predicted 
performance.  They  were  able  to  achieve  these  results  by  coupling  the  segmented  DM 
with  a  refractive  lenslet  array.  The  lenslet  improved  the  effective  fill  factor  of  the 
DM  by  focusing  the  incident  light  onto  the  centers  of  the  mirror  segments.  A  DM 
with  128  individual  elements  on  a  12  x  12  square  grid  was  used  for  this  research.  The 
combination  of  the  lenslet  and  the  segmented  DM  behaved  as  a  phase-only  modulator 
for  mirror  deflections  much  smaller  than  the  focal  length  of  the  lenslet.  [17] 

The  side  effects  associated  with  the  type  of  MEMS  DM  used  in  an  experiment 
must  be  considered  during  research.  The  number  of  actuators,  the  motion  of  resolution 
of  each  mirror  segment,  the  control  bandwidth  of  each  actuator,  and  the  actuator 
stroke  are  the  most  important  variables  that  must  be  understood  when  using  MEMS 
DMs  in  adaptive  optics  applications.  [5] 

2-4-4  Alternative  Methods. 

LC  SLMs  and  DMs  are  not  the  only  way  to  alter  a  wavefront  and  simulate 
atmospheric  turbulence.  Probably  the  most  common  method  is  a  rotating  phase 
wheel  atmospheric  turbulence  simulator.  This  device  uses  two  simple  static  phase 
plates  to  act  as  turbulence  layers.  The  phase  plates  are  created  based  on  the  concept 
of  near-index  matching.  When  they  are  rotated,  the  light  experiences  a  randomly 
evolving  index  of  refraction.  Most  of  the  important  atmospheric  parameters  can  be 
manipulated  by  changing  the  relative  positions  and  rotation  speeds  of  the  phase  plates 
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when  two  static  phase  plates  are  put  together  in  the  same  optical  system.  [30]  This  is 
not  an  ideal  method  of  turbulence  creation  because  it  takes  considerably  more  time 
to  make  simple  changes  to  the  atmosphere  than  the  other  methods  discussed.  [22] 
As  evidenced  in  any  educational  AO  laboratory,  obtaining  desired  values  for  the  at¬ 
mospheric  turbulence  parameters  requires  tedious  calculations  and  movement  of  the 
phase  wheels.  Phase  wheels  are  often  repositioned  by  hand  which  is  not  as  precise  as 
an  electronically  addressed  wavefront  control  device  such  as  a  SLM  or  DM.  SLMs  and 
DMs  are  much  more  reproducible  because  they  are  electronically  addressed.  SLMs 
also  benefit  from  being  nonmechanical  in  that  they  do  not  have  any  moving  parts 
that  require  maintenance.  Phase  wheel  systems  can  also  be  cumbersome  and  usually 
require  several  meters  of  space  on  an  optical  bench,  whereas  SLMs  and  MEMS  DMs 
are  compact  and  require  only  a  few  centimeters  of  space  in  a  laboratory. 
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III.  Theoretical  Analysis 


This  chapter  conveys  the  theoretical  approach  and  methodology  used  to  address 
the  main  ideas  of  this  research.  First,  a  theoretical  understanding  of  setting 
up  two  phase  screens  for  atmospheric  simulation  is  presented  and  discussed.  Then, 
a  detailed  analysis  of  the  minimum  atmospheric  coherence  diameter  as  a  function  of 
spatial  frequency  is  presented  and  discussed. 

3.1  Phase  Screen  Locations 

This  section  investigates  the  strongest  possible  optical  turbulence  that  can  be 
emulated  by  two  wavefront  control  devices  in  a  laboratory.  Previous  work  by  Phillips 
et  al.  [27]  concluded  that  a  broad  range  of  turbulence  parameters  could  be  represented 
using  only  two  LC  SLMs  as  phase  screens.  Two  phase  screens  are  commonly  used  in  a 
laboratory  environment  due  to  cost  limitations  and  their  ability  to  represent  a  broad 
range  of  atmospheric  conditions.  This  research  quantifies  the  maximum  strength  of 
turbulence  two  phase  screens  are  capable  of  representing. 

The  expression  for  atmospheric  coherence  diameter  is  given  in  a  discrete  form 
by  Eq.  (2.16).  When  this  equation  is  expanded  for  two  phase  screens  it  becomes 


(3.1) 


It  is  only  dependent  on  the  atmospheric  coherence  diameters  of  the  individual  phase 
screens.  The  value  r0  represents  the  distance  between  two  coherent  points  in  the 
atmosphere.  Therefore,  turbulence  is  at  a  maximum  strength  when  that  distance  is 
small.  In  the  equation  above,  and  tq2  are  independent  of  each  other  and  thus, 
the  minimum  possible  overall  atmospheric  coherence  length  occurs  when  both  phase 
screen  r0  values  are  at  a  minimum.  Assuming  both  phase  screens  are  created  with  the 
same  type  of  device,  this  limit  is  the  same  for  each  screen.  Under  these  constraints, 
Eq.  (3.1)  simplifies  to 


ro  2  /  r0iinin(Fmax), 


(3.2) 


where  Fmax  is  the  Nyquist  bandlimit  of  the  pixelated  wavefront  control  device  used 
as  a  phase  screen,  and  roimin  is  the  minimum  possible  value  for  atmospheric  coherence 
length  that  can  be  represented  with  each  phase  screen.  A  segmented  wavefront  control 
device  has  a  Nyquist  bandlimit  corresponding  to 

F  =  — 

- L  max  2J  ’ 

where  Sx  is  the  width  of  each  individual  pixel.  [13]  The  number  of  pixels  N  across  a 
device  can  be  found  with 

N=^  =  2  DFmax,  (3.4) 

Ox 

where  D  is  the  diameter  of  the  active  area  of  the  device.  Using  these  equations,  it  is 
easy  to  see  that  is  both  a  function  of  the  spatial  frequency  of  the  device  as  well 
as  the  number  of  pixels  across  the  device. 

The  other  atmospheric  parameter  of  interest  is  the  Rytov  number,  which  is 
also  known  as  the  log-amplitude  variance  under  weak  turbulence  conditions.  The 
equation  for  this  value  was  given  for  a  layered  turbulence  model  by  Eq.  (2.17)  and 
can  be  written  for  two  phase  screens  by 
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Using  the  same  rationale  as  before,  the  individual  coherence  diameters  both  become 
7Y)  .  and  the  expression  above  reduces  to 
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for  the  maximum  possible  Rytov  number.  Creating  linear  arrays  of  values  for  the 
variables  for  z\/L  and  z-zjL  leads  to  a  contour  plot  for 
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Figure  3.1:  Maximum  normalized  Rytov  number  yTxsfe  J  giyen  as  a  function  of 

different  phase  screen  locations  for  a  two-screen  setup. 


which  is  shown  in  Fig.  3.1.  This  equation  shows  that  the  maximum  Rytov  number 
is  a  function  of  r0imin,  total  propagation  distance,  wavelength,  and  also  the  location 
of  the  phase  screens  for  a  two-screen  setup.  Most  importantly,  Fig.  3.1  shows  that 
the  maximum  value  for  the  Rytov  number  can  be  found  when  both  phase  screens  are 
as  close  as  possible  to  the  beginning  of  the  propagation  distance.  This  is  possible 
because  Eq.  (3.6)  does  not  depend  on  the  total  length  of  atmosphere  each  phase 
screen  represents  Azi,  but  only  on  the  physical  position  of  each  phase  screen  Zi.  Even 
though  the  theoretical  setup  for  the  most  turbulent  atmosphere  requires  both  phase 
screens  to  be  placed  at  the  very  beginning  of  the  optical  path,  the  physical  size  of 
each  wavefront  control  device  limits  their  placement.  This  physical  limitation  may  be 


30 


insurmountable  to  achieve  the  maximum  possible  turbulence,  but  Fig.  3.1  provides 
an  estimate  for  the  maximum  Rytov  number  obtainable  given  any  laboratory  setup 
for  a  two  phase  screen  system. 

3.2  Minimum  Atmospheric  Coherence  Length 

The  previous  discussion  on  phase  screen  locations  has  revealed  that  the  mini¬ 
mum  possible  overall  value  for  r0  and  the  maximum  possible  overall  value  for  cr^  are 
functions  of  roimin.  In  order  to  truly  identify  the  maximum  turbulence  that  can  be 
created  in  a  laboratory  with  two  phase  screens,  it  is  essential  to  look  into  the  absolute 
minimum  value  for  the  atmospheric  coherence  length  of  a  phase  screen.  As  previously 
stated,  this  value  is  a  function  of  the  Nyquist  bandlimit  of  the  segmented  wavefront 
control  device  being  used.  Equation  (3.4)  shows  there  is  a  direct  relationship  between 
the  Nyquist  bandlimit  and  the  number  of  pixels  across  a  device.  Therefore,  an  inves¬ 
tigation  is  required  to  determine  the  relationship  for  r0imin  as  a  function  of  the  number 
of  pixels  across  a  device. 

The  initial  method  used  for  Ending  the  minimum  atmospheric  coherence  diam¬ 
eter  in  terms  of  spatial  frequency  is  to  begin  with  the  refractive  index  PSD,  convert 
it  to  overall  phase  variance,  and  then  solve  that  for  r0.  The  first  PSD  investigated  is 
the  Tatarskii  spectrum  as  given  by 


where  =  5.92 /l0.  [1]  This  spectrum  is  used  because  it  includes  the  inner  scale 
lo  which  limits  the  PSD  at  higher  spatial  frequencies.  When  written  in  terms  of 
atmospheric  coherence  diameter,  Noll  [25]  shows  how  this  equation  reduces  to 
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where  v  is  the  radial  polar  component  of  the  three-dimensional  linear  spatial  frequency 
vector  v  and  is  measured  in  waves/meter.  The  phase  PSD  is  not  the  most  convenient 
parameter  to  evaluate  so  the  overall  phase  variance  associated  with  the  wavefront  is 
used  instead.  Phase  variance  can  be  found  by  integrating  the  PSD  and  is  given  by 


POO 

<otai  =  27T  ®(v)vdv.  (3.10) 

Jo 

This  equation  gives  the  phase  variance  (measured  in  rad2)  for  all  spatial  frequen¬ 
cies  from  zero  to  infinity.  This  research,  however,  is  interested  in  identifying  how 
limitations  in  spatial  frequency  effect  the  range  of  turbulence  parameters.  In  order 
to  represent  this  constraint,  Eq.  (3.10)  can  only  be  evaluated  from  Fmax  to  infinity. 
Constraining  the  limits  of  integration  allows  one  to  find  only  the  phase  variance  not 
simulated  by  the  segmented  wavefront  control  device.  The  overall  range  of  repre¬ 
sentable  phase  variance  by  the  device  can  be  found  with 
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In  order  to  evaluate  r0  in  realistic  terms  and  remove  extraneous  variables,  a 
level  of  fidelity  must  be  established.  The  level  of  fidelity  can  best  be  created  by 
setting  an  error  threshold  of  missing  phase  variance  that  is  acceptable  as  a  result  of 
using  a  segmented  wavefront  control  device.  This  error  threshold  is  a  level  of  overall 
missing  phase  standard  deviation  cx0misBing  (measured  in  waves)  that  is  acceptable  to 
the  researcher  and  is  represented  by 
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When  measured  in  percentage  of  phase  standard  deviation  missing,  this  level  isl  — 
(%present)(o>totai);  where  %present  is  the  threshold  level  determined  by  the  percentage  of 
the  overall  phase  standard  deviation  that  the  wavefront  control  device  must  be  capable 
of  emulating.  This  research  investigates  the  threshold  of  percent  phase  standard 
deviation  present  at  the  83%,  90%,  95%,  96%,  97%,  98%,  and  99%  levels.  This  range 
should  cover  the  needs  of  most  researchers  while  providing  a  detailed  outlook  on  how 
the  varying  requirements  of  researchers  alter  the  numbers  and  types  of  devices  that 
will  satisfy  those  needs.  The  terms  “phase  error”  and  “missing  phase”  will  be  used 
interchangeably  throughout  this  investigation. 


3.2.1  Unit  Analysis. 

Now  that  a  threshold  has  been  established,  it  is  instructive  to  perform  a  unit 
analysis  of  the  equations  that  are  to  be  manipulated.  First,  it  will  be  verified  that 
Eq.  (3.10)  actually  yields  phase  variance  in  rad2.  This  analysis  begins  with  the  phase 
structure  function  D^R)  which  is  given  by  Noll  [25]  to  be 

/  D  \  5/3 

D^R)  =  6.88  (  — J  ,  (3.15) 

and  is  measured  in  rad2.  This  function  is  related  to  the  Wiener- Khintchine  theorem 
given  in  Eq.  (2.6)  by  [25] 


D*(R)  =  2  /f 
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This  equation  can  be  used  to  determine  the  units  of  the  phase  PSD  <h(z/).  When 
represented  as  units,  Eq.  (3.16)  becomes 


,9  t  .  /^/  xn  r„  /  /  m  /  waves  v 

rad  =  Umts(<3>(z/))  [1  —  cos(m/m)]  - J  , 


33 


where  the  units  in  the  cosine  divide  out.  This  results  in  the  units  of  the  PSD  given 


by 
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A  simple  substitution  of  this  expression  into  Eq.  (3.10)  yields 
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Thus  confirming  that  Eq.  (3.10)  yields  phase  variance  in  rad2. 

In  order  to  compare  the  phase  variance  to  the  threshold  value  for  phase  standard 
deviation,  these  units  are  converted  to  waves.  The  conversion  from  units  of  rad2  to 
waves  is  given  by 

Vrad2 

waves  =  - . 

2tt 

The  units  of  these  equations  must  be  carefully  tracked  so  there  is  no  confusion  in  the 
final  results.  For  this  research,  phase  variance  is  given  in  rad2  and  phase  standard 
deviation  is  always  measured  in  waves. 


3.2.2  Evaluation  Over  a  Single  Point  in  Space. 

With  the  unit  verification  complete,  it  is  possible  to  move  forward  and  evaluate 
the  atmospheric  coherence  diameter  as  a  function  of  the  number  of  pixels  across  a 
wavefront  control  device.  Equations  (3.9)  and  (3.14)  can  be  combined  to  form 
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A  simple  change  of  variables  must  be  performed  to  reduce  this  equation  to 
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where  /  =  (/0/0.942)2i/2  and  /  (Fmax)  =  (/0/0.942)2F^ax.  The  integral  portion  of  this 
equation  can  be  recognized  as  an  upper  incomplete  gamma  function  which  further 
reduces  the  equation  to 

4™L„  =  O.fKSr^3  (-g, /(*■»»,))  ■  (3.19) 


An  incomplete  gamma  function  cannot  be  evaluated  with  a  negative  value  in  its 
argument,  so  this  was  transformed  using  a  simple  identity  from  Arfken  and  Weber  [2] 
given  by 

T(a  —  l,  x)  =  aT(a,  x)  +  xae~x .  (3.20) 


Using  this  identity  to  evaluate  Eq.  (3.19)  leads  to 
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When  the  expression  for  /  (Fmax)  is  substituted  back  into  this  equation  and  it  is  solved 
for  r0,  it  becomes 
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(3.22) 

This  equation  provides  the  sought  after  expression  for  the  minimum  atmospheric 
coherence  length  in  terms  of  the  Nyquist  bandlimit  of  a  segmented  wavefront  control 
device.  It  can  be  manipulated  to  express  roimin  as  a  function  of  the  number  of  pixels 
across  a  device  using  the  relationship  between  N  and  the  Nyquist  bandlimit  given  by 
Eq.  (3.4).  When  realistic  values  are  chosen  for  the  phase  variance  error,  inner  scale, 
and  overall  aperture  diameter,  Eq.  (3.22)  can  be  plotted  as  shown  in  Fig.  3.2.  This 
figure  shows  that  Eq.  (3.22)  provides  the  minimum  possible  atmospheric  coherence 
diameter  as  the  Nyquist  bandlimit  of  a  device  increases. 

Unfortunately  there  is  no  obvious  practical  application  of  Eq.  (3.22)  at  this 
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Figure  3.2:  Minimum  possible  ro  value  as  a  function  of  the  Nyquist  bandlimit  of  a 

wavefront  control  device  for  D  =  10  mm,  Iq  =  45/im,  and  <70error  =  1/10  waves. 

point  in  time.  The  amount  of  phase  standard  deviation  missing  as  a  result  of  using  a 
pixelated  device  is  always  given  as  the  total  phase  standard  deviation  error  over  the 
entire  active  aperture  of  the  device.  Regrettably,  Eq.  (3.22)  only  gives  the  relationship 
for  a  single  point  on  that  aperture.  In  order  to  make  the  value  obtained  from  this 
equation  useful,  a  complicated  correlation  matrix  would  have  to  be  created  for  the 
phase  at  every  point  on  the  aperture.  Currently,  methods  do  not  allow  the  creation 
of  such  an  autocorrelation  matrix.  However,  future  research  is  anticipated  to  provide 
the  necessary  matrix.  The  analytic  expression  in  Eq.  (3.22)  cannot  be  verified  until 
future  research  comes  to  fruition.  Therefore,  another  method  must  be  found  for  the 
relationship  between  r0imin  and  FmXY  which  is  investigated  in  Section  3.2.3. 
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3.2.3  Evaluation  Over  an  Aperture. 

The  search  for  a  method  to  evaluate  the  total  phase  variance  over  an  aperture 
of  diameter  D  ends  with  the  use  of  Zernike  modes.  As  stated  in  Section  2.3,  a  Zernike 
polynomial  can  be  written  as  a  function  of  the  active  aperture  of  the  wavefront  control 
device.  Also,  the  total  phase  of  a  wavefront  can  be  represented  as  a  weighted  sum  of 
Zernike  modes.  Similarly,  the  total  phase  variance  of  a  wavefront  over  an  aperture 
can  be  represented  as  a  sum  of  the  phase  variance  associated  with  each  Zernike  mode 
over  the  aperture,  as  demonstrated  in  Eq.  (2.23).  [32] 

The  overall  missing  phase  variance  of  a  segmented  wavefront  control  device 
stems  from  its  limitations  on  two  fronts.  First,  the  Nyquist  bandlimit  of  the  device 
limits  the  overall  phase  it  can  represent.  Second,  in  order  to  analytically  represent 
missing  phase  variance  as  a  function  of  spatial  frequency,  the  number  of  Zernike 
polynomials  that  can  be  added  to  a  device  must  also  be  limited.  Zernike  modes  can 
reproduce  true  turbulence  when  they  are  summed  from  the  first  mode  (piston)  to 
infinity.  [25]  However,  the  best  way  to  reproduce  turbulence  as  a  function  of  spatial 
frequency  is  to  write  it  as  a  finite  sum  of  Zernike  modes.  The  modes  excluded  from 
this  sum  represent  the  part  of  the  total  missing  phase  variance  from  finite  modes. 
An  analytical  expression  for  r o,.miTi  as  a  function  of  Nyquist  bandlimit  can  not  derived 
for  the  Zernike  mode  approach.  A  numerical  approach  works,  but  only  for  a  finite 
number  of  Zernike  modes.  This  research  only  focuses  on  the  higher  order  modes  that 
a  device  can  emulate.  Piston,  tip,  and  tilt  can  all  be  easily  represented  in  turbulence 
with  the  use  of  a  fast  steering  mirror  (FSM).  Tip  and  tilt  represent  approximately 
87%  of  the  total  phase  variance  in  a  turbulent  path.  [32]  This  causes  much  of  the 
throw  of  the  actuators  on  a  complex  device  such  as  a  MEMS  DM  to  be  wasted  on 
producing  the  first  three  modes  and  limits  the  total  range  available  to  represent  the 
higher  order  modes  where  i  is  equal  to  four  or  greater.  For  this  reason,  only  the  phase 
variance  associated  with  higher  order  modes  has  been  investigated  and  it  is  assumed 
that  all  of  the  phase  from  piston,  tip,  and  tilt  are  represented  with  another  device 
such  as  a  FSM. 
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3.2.3. 1  Fidelity  Levels. 


A  level  of  fidelity  must  be  chosen  to  determine  the  amount  of  overall  phase 
variance  missing  from  the  combined  effects  of  pixelation  and  finite  modes.  This  level 
represents  an  acceptable  limit  to  fidelity.  Once  this  total  fidelity  level  is  chosen, 
it  is  budgeted  between  the  amount  missing  due  to  finite  modes  and  the  amount 
missing  from  the  Nyquist  bandlimit  of  the  device.  Since  the  amount  of  phase  standard 
deviation  missing  from  pixelation  is  dependent  on  the  number  of  modes,  the  level 
from  finite  modes  must  be  chosen  first.  For  example,  if  the  researcher  determines 
that  it  is  acceptable  to  represent  only  95%  of  the  total  phase  standard  deviation, 
the  amount  present  by  having  finite  modes  will  be  97.5%  of  the  total.  Then  the 
amount  representable  by  using  a  pixelated  device  will  be  97.5%  of  the  phase  standard 
deviation  left  after  losing  some  due  to  finite  modes.  These  values  are  given  by 

%  of  Total  Present  =  (%  of  Total  Present  due  to  Finite  Modes) 

x  (%  of  Modes  Present  due  to  Pixelation) .  (3.23) 

In  an  effort  to  simplify  these  two  fidelity  levels,  the  percentage  of  the  total  phase 
standard  deviation  present  due  to  finite  modes  was  set  equal  to  the  percentage  present 
after  pixelating  the  total  mode  phase  standard  deviation.  This  reduces  the  above 
equation  to 


%  of  Total  Present  =  (%  Present  due  to  Modes  or  Pixelation)2  .  (3.24) 


An  example  is  given  to  clarify  this  confusing  concept.  According  to  Sasiela  [32], 
the  total  amount  of  phase  standard  deviation  that  exists  in  a  wavefront  after  piston, 
tip,  and  tilt  have  been  removed  (PTR)  is  given  by 


(3.25) 


38 


as  measured  in  waves.  When  D/r0  is  normalized  out  of  this  equation,  the  total  PTR 
phase  standard  deviation  cr0PTR  of  a  wavefront  is  approximately  0.058260  waves.  If  the 
researcher  can  accept  a  total  loss  of  5%  of  phase  standard  deviation  over  an  aperture, 
that  leaves  (0.95)  (0.058260)  ~  0.055347  waves  of  phase  standard  deviation  that  must 
be  represented  with  the  wavefront  control  device.  The  percentage  present  after  only 
a  finite  number  of  modes  has  been  represented  can  be  found  using  Eq.  (3.24)  and  is 
given  as  (0.95)1/2  ~  0.975.  This  leads  to  (0.975) (0.058260)  ~  0.056804  waves  of  total 
phase  standard  deviation  present  after  summing  to  a  finite  number  Zernike  modes. 
The  total  phase  standard  deviation  of  that  value  present  after  considering  the  Nyquist 
bandlimit  of  the  device  is  also  97.5%  as  given  by  the  fidelity  level  equations.  This 
leads  to  (0.975)  (0.056804)  ~  0.0553839  waves  of  total  PTR  phase  standard  deviation 
that  a  device  can  represent  after  losing  some  to  finite  modes  and  some  to  the  spatial 
frequency  of  the  device.  This  value  is  approximately  95%  of  the  total  PTR  phase 
standard  deviation  of  a  wavefront  as  given  by  Sasiela  [32]  in  Eq.  (3.25).  Therefore, 
this  example  shows  how  levels  of  fidelity  are  evaluated  in  this  research.  A  total  loss 
of  5%  can  be  found  by  losing  2.5%  of  the  total  phase  standard  deviation  present  due 
to  finite  modes  and  then  losing  another  2.5%  of  that  due  to  the  Nyquist  bandlimit  of 
a  wavefront  control  device. 

This  example  and  analysis  method  leads  to  the  threshold  levels  shown  in  Ta¬ 
ble  3.1  for  both  finite  modes  and  pixelation  of  a  device,  given  a  threshold  for  the 
overall  normalized  phase  standard  deviation  present  in  a  device.  Several  common 
threshold  levels  are  shown  in  this  table  as  well  as  the  83%  threshold  level  which  is 
also  in  all  subsequent  tables  in  this  research  to  support  an  example  in  Chapter  IV. 

3. 2. 3. 2  Phase  Standard  Deviation  from  Finite  Zernike  Modes. 

It  has  been  shown  in  Section  3.2.3. 1  that  the  total  phase  standard  deviation 
missing  depends  on  the  total  number  of  Zernike  modes  needed  to  reach  a  certain 
fidelity  level  of  total  phase  standard  deviation  present  after  removing  piston,  tip,  and 
tilt.  To  find  this,  an  equation  is  needed  for  the  total  phase  variance  present  in  a  given 
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Tabic  3.1:  Threshold  percentages  for  finite  modes  and  pixelation  of  a  device  given  a 

threshold  for  the  overall  normalized  phase  standard  deviation  present  in  a  device. 


Percentage  of  total 
a<t>PTR  Present 

Percentage  present  from 
finite  modes  %presmodeB 

Percentage  present  from 
pixelation  %presDixels 

99% 

99.5% 

99.5% 

98% 

99.0% 

99.0% 

97% 

98.5% 

98.5% 

96% 

98.0% 

98.0% 

95% 

97.5% 

97.5% 

90% 

94.9% 

94.9% 

83% 

91.1% 

91.1% 

Zernike  mode  cr^..  Fortunately,  Roggemann  and  Welsh  [31]  offer  this  expression  as 


4.5  x  10~4tt8/3 


r(ir)r(<3-ir) 

R(¥)]2r  (0  +  $)’ 


(3.26) 


where  a2  is  the  variance  of  the  Zernike  coefficients,  i  is  the  Zernike  mode  in  question, 
/3  =  n  +  1,  and  a  =  16/3.  Recall  from  Section  2.3  that  n  and  m  are  directly  correlated 
to  the  ith  Zernike  mode,  as  referenced  in  Table  2.1.  Equation  (3.26)  does  not  allow  for 
the  next  step  in  the  process.  Therefore,  an  equation  is  needed  that  provides  the  phase 
variance  present  in  a  Zernike  mode  as  a  function  of  spatial  frequency.  Sasiela  [32]  gives 
such  an  equation  with 

uj.  =  0.2073 k2  f  dzC2n(z) 

Jo 

where  P{ 7,  k,  z)  specifies  the  type  of  optical  source  and  is  zero  due  to  no  actual  prop¬ 
agation  taking  place.  Further,  f(n)  =  ac-11'3,  as  given  by  the  Kolmogorov  spectrum, 
and  F( 7,  k)  is  a  filter  function  used  to  represent  the  individual  Zernike  modes  given 


I  dKf(K)  cos2  [P(7,  «,^)]F(7«),  (3.27) 
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by 


/ 


2  cos  2(m0) 


F(7k)  =  Fm,n(K)„ 


1 


where  m  and  n  are  the  azimuthal  and  radial  orders,  respectively,  and  Jn+1  is  a  Bessel 
function  of  the  first  kind.  This  analysis  requires  the  use  of  the  Kolmogorov  spectrum 
and  discontinues  a  PSD  that  is  dependent  on  the  inner  scale  because  an  analytic 
solution  cannot  be  found  for  the  total  phase  variance  in  a  given  mode  with  these  more 
complex  equations.  Equation  (3.27)  can  be  written  in  terms  of  r0  after  substituting 
Eq.  (2.15).  Also,  when  all  modes  are  being  considered  Fm^n(hi)x  and  Fm>n{K)y  can  be 
combined  so  that  the  trigonometric  functions  drop  out  due  to  the  identity 


sin2  (a?)  +  cos2  (a;)  =  1. 


(3.29) 


The  result  is  a  simpler  equation  for  the  total  phase  variance  present  in  each  mode  as 
a  function  of  r0  and  spatial  frequency  given  by 


(3.30) 


where  once  again  j3  =  n  +  1  and  a  =  16/?. 

Before  investigating  how  this  equation  can  lead  to  a  relationship  between  the 
minimum  atmospheric  coherence  diameter  and  the  Nyquist  bandlimit  of  a  device,  it 
is  important  to  verify  Eq  (3.30).  If  a  device  is  continuous,  the  limits  of  integration 
shown  in  Eq.  (3.30)  are  correct  and  they  represent  the  total  phase  variance  present  in 
the  ith  Zernike  mode.  This  value  should  be  equal  to  that  given  by  Eq.  (3.26)  which 
implicitly  includes  all  spatial  frequencies.  In  order  to  compare  these  two  equations,  an 
analytic  form  must  be  found  for  Eq.  (3.30).  The  help  of  Mellin  transforms  A 4{f(t)]  s} 
can  be  used  to  obtain  such  a  solution.  Poularikas  [28]  gives  the  scaling  property  of  a 
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Mellin  transform  as 


M{f(at);s} 


f(at)ts  ldt  =  a 


poo 

~s  f  f(x)xs-1dx  =  a~sF(s ) 
o 


(3.31) 


Sasiela  [32]  makes  this  property  useful  with  the  Mellin  transform  pair  for  the  square 
of  a  Bessel  function  of  the  first  kind  which  is  given  by 


i  r(f  +  /?)r(i-f) 
2v/irp+i-i)r(i-f) 


(3.32) 


Using  the  relationships  expressed  in  Eqs.  (3.31)  and  (3.32)  to  simplify  Eq.  (3.30)  leads 
to 

<  =  °-49007(2r  0 5/3  r(|}r(/3  + j)’  (3'33) 

which  is  another  analytic  expression  for  the  total  phase  variance  present  in  the  ith 
Zernike  mode.  The  difference  is  that  Eq.  (3.33)  is  derived  from  an  equation  given 
as  a  function  of  spatial  frequency.  This  expression  looks  similar  to  that  given  by 
Roggemann  and  Welsh  in  Eq.  (3.26).  However,  to  verify  that  both  expressions  yield 
the  same  values,  they  can  be  normalized  so  that  D /r$  is  removed  from  the  equations 
and  plotted  on  the  same  graph.  This  graph  is  shown  in  Fig.  3.3  and  it  is  visually 
apparent  that  Eqs.  (3.26)  and  (3.33)  yield  the  same  values.  Now  that  a  method 
for  expressing  the  total  phase  standard  deviation  as  a  function  of  modes  has  been 
presented,  the  number  of  modes  corresponding  to  the  thresholds  defined  in  Table  3.1 
can  be  found.  These  modes  are  shown  in  Fig.  3.4  and  Table  3.2. 


3. 2. 3. 3  Phase  Standard  Deviation  from  Finite  Number  of  Pixels. 

An  expression  for  phase  standard  deviation  as  a  function  of  spatial  frequency 
can  be  investigated  with  the  number  of  modes  known  for  each  fidelity  threshold  of 
interest  to  this  research.  The  verification  of  Eq.  (3.33)  implies  that  the  equation  it 
was  derived  from  is  also  correct.  Equation  (3.30)  must  be  normalized  so  that  D/r0 
can  be  removed  from  the  integral.  This  is  accomplished  by  converting  the  spatial 
frequency  variable  k  to  a  dummy  variable  for  normalized  linear  spatial  frequency  r), 
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Figure  3.3:  Comparison  of  total  normalized  phase  standard  deviation  as  a  function  of 

modes  equations  from  Roggemann  and  Welsh  and  this  research’s  derivation 


where 


Tj  =  k- 


D 


(3.34) 


This  substitution  allows  the  phase  variance  present  in  each  mode  to  be  given  by 

(3.35) 


o. 


5/3  /*oc 

tqJ  Jo 


2I  =  0.98014(2)~11/37ra  (  — ]  /  [Jp  (r/)]2  r]~u/3dr]: 


which  is  only  a  function  of  the  variable  ij  and  is  directly  related  to  the  total  number 
of  pixels  on  a  device  when  D/tq  is  normalized  out  of  the  equation.  The  expression 


^max  =  (3.36) 

converts  the  total  number  of  pixels  on  a  device  to  the  maximum  value  of  the  variable 
r)  to  determine  how  much  phase  variance  is  present  on  a  device  with  a  limited  number 
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Figure  3.4:  Number  of  modes  required  to  reach  a  variety  of  thresholds  of  total  phase 

standard  deviation  present  after  using  a  discrete  number  of  Zernike  modes. 


of  pixels.  Substituting  this  into  the  upper  limit  of  integration  in  Eq.  (3.35)  gives 

/  T~\  \  5/3  /*?7niax 

°ii(v)  =  0.98014(2)-11/37ra  (— J  J [Jp  (rj)]2  rj~u/3drj,  (3.37) 

which  is  an  expression  for  the  total  amount  of  phase  variance  present  in  a  given  mode 
due  to  the  finite  pixels  on  a  wavefront  control  device. 

Finally,  the  amount  of  phase  variance  present  on  a  segmented  wavefront  control 
device  with  a  limited  number  of  pixels  after  a  finite  number  of  higher  order  Zernike 
modes  have  been  summed  can  be  expressed  by 

total  modes  present 

^pixels  (V)  =  Yl  (3-38) 

i— 4 
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Tabic  3.2:  Number  of  modes  required  to  reach  certain  thresholds  of  total  phase  standard 

deviation  present  after  using  a  discrete  number  of  Zernike  modes. 


Percentage  of  total 
a<PpTR  present 

Percentage  present 
from  finite  modes 

Number  of  modes 
required  to  reach  threshold 

99% 

99.5% 

457 

98% 

99.0% 

218 

97% 

98.5% 

139 

96% 

98.0% 

101 

95% 

97.5% 

78 

90% 

94.9% 

35 

83% 

91.1% 

20 

Earlier  in  this  section,  it  was  stated  that  piston,  tip,  and  tilt  are  not  included  in 
this  sum  because  they  are  easy  to  include  with  another  device  and  interrupt  the  true 
number  of  higher  order  modes  that  a  complex  wavefront  control  device  can  produce. 

Figure  3.5  shows  the  plot  of  Eq.  (3.38).  This  plot  shows  the  total  normalized 
phase  standard  deviation  present  as  a  function  of  spatial  frequency  for  several  levels 
of  fidelity  after  limiting  the  number  of  Zernike  modes  present  on  a  wavefront.  Since 
different  threshold  levels  require  a  different  number  of  Zernike  modes  to  be  summed  in 
Eq.  (3.38),  the  plots  are  slightly  different  for  each  threshold  level.  The  steepest  slope 
is  seen  at  the  lowest  threshold  level,  where  there  are  the  fewest  Zernike  polynomials 
summed.  This  is  intuitive  because  fewer  modes  mean  a  lower  total  phase  standard 
deviation  present,  and  therefore,  it  takes  fewer  77  samples  to  sum  to  that  number. 
The  number  of  pixels  required  for  a  wavefront  control  device  for  each  threshold  level 
shown  in  Table  3.1  are  given  in  Table  3.3.  Figures  3.3,  3.4,  and  3.5  are  the  result  of 
the  Matlab®  code  shown  in  Appendix  A. 

The  results  shown  in  Table  3.3  provide  the  total  number  of  modes  and  number 
of  pixels  a  device  must  have  to  successfully  create  a  certain  percentage  of  the  tur¬ 
bulence  associated  with  all  higher  order  Zernike  modes.  It  suggests  that  to  recreate 
the  atmosphere  from  all  higher  order  modes  with  only  a  five  percent  loss,  a  wavefront 
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Figure  3.5:  Number  of  pixels  required  to  reach  a  variety  of  threshold  levels  of  total  phase 
standard  deviation  present  after  using  a  Nyquist  bandlimited  device  and  only  summing  over 
a  discrete  number  of  Zernike  modes. 


control  device  with  only  six  or  seven  pixels  across  would  be  necessary.  This  value 
seems  extremely  low.  Thus  it  is  important  to  verify  these  results  with  simulation. 


3. 2. 3. 4  Minimum  Theoretical  Atmospheric  Coherence  Length. 

Simulation  verification  is  shown  in  Chapter  IV;  however,  it  is  first  necessary 
to  complete  the  analysis  of  finding  the  minimum  possible  value  for  ro  that  can  be 
represented  by  all  higher  order  modes  with  a  specific  wavefront  control  device.  This 
expression  comes  from  using  Eq.  (3.25)  for  a  specific  device.  The  minimum  r0  is 
the  value  at  which  the  achieved  phase  standard  deviation  is  close  enough  to  the 
desired  phase  standard  deviation.  Close  enough  is  defined  by  the  desired  fidelity 
limits  previously  described  in  Section  3.2.3. 1.  This  value  should  not  be  confused  with 
the  effective  r0  of  the  turbulence  on  the  device.  Standard  MEMS  DMs  have  five  to 
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Table  3.3:  Number  of  pixels  theoretically  required  to  reach  certain  threshold  levels  of 

total  phase  standard  deviation  present  after  using  a  discrete  number  of  Zernike  modes. 


Percentage  of  total 
a<t>PTR  present 

Finite  mode 

%threshold 

#  of  modes 
required 

Finite  pixel 

%threshold 

7max 

#  of  pixels 
required 

99% 

99.5% 

457 

99.5% 

24 

15.28 

98% 

99.0% 

218 

99.0% 

16 

10.19 

97% 

98.5% 

139 

98.5% 

13 

8.28 

96% 

98.0% 

101 

98.0% 

11 

7.00 

95% 

97.5% 

78 

97.5% 

10 

6.37 

90% 

94.9% 

35 

94.9% 

7 

4.46 

83% 

91.1% 

20 

91.1% 

6 

3.82 

twelve  mirror  segments  across  the  device,  so  an  example  is  presented  using  a  seven 
pixel  device  to  better  understand  this  concept. 

A  7  x  7  pixel  device  is  able  to  represent  96%  of  the  total  PTR  phase  standard 
deviation  according  to  Table  3.3.  Therefore,  four  percent  of  the  total  is  missing  when 
a  device  with  only  seven  pixels  across  is  used.  The  minimum  possible  value  for  r$  that 
this  device  can  represent  is  a  function  of  how  many  waves  of  phase  standard  deviation 
a  researcher  is  willing  to  accept  as  missing.  This  expression  is  given  by 


C%PTRt, 


%missing 


2vr  ) 


(3.39) 


where  is  the  number  of  waves  of  missing  PTR  phase  standard  deviation 

1  ^missing 

that  a  researcher  can  accept  and  %missmg  is  the  percentage  of  missing  PTR  phase 
standard  deviation  as  a  result  of  the  number  of  pixels  on  a  device.  Solving  Eq.  (3.39) 
for  r0  .  yields 

umin  «/ 


r0 


=  D 


^0.134 

%  missing  \ 

2t r 

k  ^PTR  •  •  / 

\  T±  -*■  A  ‘'missing  / 

6/5 


(3.40) 


Equation  (3.40)  indirectly  relates  the  minimum  atmospheric  coherence  length  to  the 
number  of  pixels  across  a  device  through  the  use  of  Eqs.  (3.36),  (3.37),  and  (3.38). 
These  equations  all  lead  to  a  value  for  %rTlissi-ng  for  a  specific  number  of  pixels  across 
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[waves] 


Figure  3.6:  Minimum  atmospheric  coherence  length  as  a  function  of  missing  PTR  phase 

standard  deviation  for  theoretical  results  for  a  7  x  7  pixel  device  with  an  aperture  diameter 
D  =  10  mm. 


a  device  which,  in  turn,  leads  to  the  minimum  possible  r0min  value  for  that  device.  In 
the  example,  N  =  7  and  %missing  =  4%.  Equation  (3.40)  is  shown  as  a  function  of 
ctj,  in  Fig.  3.6  when  the  aperture  diameter  is  set  arbitrarily  as  D  =  10  mm. 

'Tt'  i  rvmissing 

Figure  3.6  shows  that,  in  theory,  a  7  x  7  pixel  device  with  an  aperture  of  D  =  10  mm 
can  represent  a  minimum  value  for  ro  =  1.736  mm  when  the  researcher  is  willing  to 
have  1/100  waves  of  missing  PTR  phase  standard  deviation.  An  even  smaller  value 
for  r0min  can  be  obtained  if  the  researcher  can  accept  a  larger  amount  of  missing 
PTR  phase  standard  deviation.  These  results  show  that  the  maximum  strength  of 
turbulence  that  a  phase  screen  can  represent  is  dependent  on  the  requirements  of  the 
researcher.  The  results  gained  from  this  analysis  will  be  discussed  in  Chapter  V. 
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IV.  Simulation  and  Analysis 


This  chapter  discusses  how  the  results  from  the  previous  chapter  were  verified 
with  computer  simulation,  ft  begins  with  the  creation  and  pixelation  of  phase 
screens  and  then  moves  on  to  how  phase  variance  is  extracted  from  a  computer  model. 
Finally,  the  results  from  the  simulation  are  presented. 

4-1  Creating  Phase  Screens 

The  first  step  in  simulating  the  theoretical  results  obtained  in  the  previous 
chapter  is  to  create  atmospheric  phase  screens.  The  goal  of  these  simulations  is  to 
verify  the  number  of  pixels  across  a  device  required  to  represent  a  certain  threshold 
level  of  phase  standard  deviation  present  in  the  emulated  turbulence  over  an  aperture. 
This  goal  is  limited  to  the  characteristics  of  a  phase  screen  on  a  single  wavefront 
control  device.  This  eliminates  all  need  for  simulating  optical  propagation  through 
the  atmosphere.  All  that  is  necessary  in  this  verification  is  to  create  a  random  number 
of  phase  screens,  compute  the  statistics  of  those  screens,  and  compare  those  to  the 
theoretical  values  calculated  in  the  previous  chapter. 

One  of  the  most  accepted  ways  to  create  a  turbulent  phase  screen  in  simulation 
is  to  use  the  Fourier  series  method  discussed  in  Section  2.2.2.  This  method  allows 
for  near-continuous  phase  screens  that  incorporate  a  large  number  of  Fourier  series 
coefficients,  which  produces  more  realistic  atmospheric  phase  screens  than  is  possible 
by  summing  a  finite  number  of  Zernike  modes.  The  Fourier  series  method  produces 
phase  screens  in  simulation  that  are  similar  to  recreating  phase  screens  by  summing 
over  a  near-infinite  number  of  Zernike  modes. 

This  research  capitalized  on  previous  work  performed  by  Dr.  Eric  Magee  of  MZA 
Associates,  Inc.  and  used  several  of  his  Matlab®  functions  to  create  the  necessary 
phase  screens.  Appendix  B  contains  the  Matlab®  hie  written  to  create  phase  screens 
via  the  Fourier  series  method.  This  listing  creates  250  random  phase  screens  over 
an  aperture  of  unit  radius  all  with  values  for  D/r§  =  5.  It  uses  the  Kolmogorov 
turbulence  spectrum  discussed  in  Section  2.1  to  produce  screens  with  128  Fourier 
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Figure  4.1:  Random  individual  realization  of  a  piston/tilt  removed  phase  screen  created 

using  the  Fourier  series  method. 

series  coefficients.  A  large  number  of  realizations  is  needed  to  accurately  calculate 
the  statistics  of  the  phase  screens.  Two  hundred  and  fifty  random  phase  screen 
realizations  offer  a  large  enough  sample  base  to  provide  good  statistical  values  of  the 
phase  variance.  D/r0  can  be  set  to  any  number  because  it  is  eventually  normalized 
out  of  the  simulation  results.  Enough  Fourier  series  coefficients  are  needed  to  create 
realistic  turbulence  in  the  phase  screens.  The  value  chosen  allows  for  realistic  near- 
continuous  phase  screens  without  taking  too  much  computational  effort.  First,  this 
listing  calls  the  Matlab®  function  mkPSD.thesis  .  m  which  creates  the  individual  power 
spectral  density  realizations  for  the  phase  screens  to  be  created  with  the  Fourier 
series  method.  Next,  Listing  B.l  calls  the  Matlab®  function  FSCoeff  .thesis  .  m  which 
calculates  the  Fourier  series  coefficients  from  the  individual  PSD  realizations.  After 
that,  the  listing  calls  the  Matlab®  function  FSscreenReal.thesis  .m  which  actually 
generates  the  individual  phase  screen  realizations.  The  three  Matlab®  functions  called 
were  all  provided  by  Dr.  Magee  and  may  be  available  upon  request  via  email  at 
emagee@mza.com. 
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After  an  individual  phase  screen  has  been  created,  Listing  B.l  calculates  the 
Zernike  coefficients  for  piston,  tip,  and  tilt.  This  allows  the  program  to  then  remove 
the  first  three  Zernike  modes  from  the  phase  screen  realization.  It  was  shown  in 
Section  2.3  that  Zernike  polynomials  can  be  removed  from  a  phase  screen  through 
simple  subtraction  due  to  each  mode’s  orthogonality  over  a  circular  aperture.  A  PTR 
phase  screen  is  necessary  because  this  research  focuses  on  the  ability  of  a  wavefront 
control  device  to  create  atmosphere  with  only  the  higher  order  modes  included.  A 
FSM  can  easily  reproduce  the  effects  of  the  first  three  Zernike  modes  as  previously 
discussed  in  Section  3.2.3.  The  final  result  of  this  Matlab®  code  is  250  random 
realizations  of  piston/tilt  removed  phase  screens.  One  of  these  realizations  is  shown 
in  Fig.  4.1.  This  figure  is  representative  of  all  250  realizations  since  each  one  is  created 
randomly. 

f.2  Pixelation  of  Phase  Screens 

The  next  step  in  the  simulation  process  is  to  down-sample  each  phase  screen. 
Limitations  in  Matlab®’s  processing  capabilities  limit  the  size  of  each  phase  screen 
to  512  x  512  pixels.  This  limitation  means  that  the  near  continuous  sample  shown  in 
Fig.  4.1  is  actually  a  512  x  512  pixelated  phase  screen.  In  order  to  find  the  effect  of 
pixelation  on  the  standard  deviation  across  many  realizations  of  a  phase  screen,  it  is 
necessary  to  create  an  array  of  pixelated  phase  screens  for  which  to  find  the  statistics. 

The  Matlab®  code  that  down-samples  the  phase  screen  realizations  can  be 
viewed  in  Appendix  C.  It  is  essentially  an  averaging  function  that  down-samples 
the  phase  screen  by 

jv»=.  =  Ay  (4,i) 

n 

where  Anew  is  the  number  of  pixels  desired  across  the  down-sampled  phase  screen, 
Norig  is  the  number  of  pixels  across  the  original  phase  screen,  and  n  is  the  value  that 
satisfies  this  equation.  The  Matlab®  code  maps  the  original  phase  screen  onto  a 
AneW  x  Anew  grid  with  each  square  inside  the  grid  composed  of  n  x  n  pixels  of  the 
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Figure  4.2:  fVnew  x  Nnew  =  4x4  grid  showing  each  square  of  n  x  n  =  16  pixels  that  are 
averaged  to  down-sample  a  Nor\g  x  Aorig  =  16  x  16  pixel  phase  screen. 

original  phase  screen.  This  grid  is  shown  in  Fig.  4.2  for  an  example  where  iVorig  =  16 
is  down-sampled  to  Nnew  =  4  with  anxn  =  4x4  grid.  Next,  each  pixel  in  the  n  x  n 
square  is  replaced  with  the  average  of  all  the  pixels  inside  the  square.  Repeating  this 
method  at  every  square  on  the  grid  leads  to  a  phase  screen  that  is  actually  Nmig  pixels 
across,  but  has  been  down-sampled  to  simulate  a  screen  that  is  only  iVnew  pixels  across. 
Each  down-sampled  screen  remains  Nov\g  x  Nm ;g  pixels  which  enables  it  to  better  map 
to  the  matrix  sizes  of  the  other  variables  in  the  code.  Also,  keeping  the  pixelated 
screens  the  original  size  enables  them  to  easily  be  placed  on  a  SLM  in  the  future  if 
experimental  verification  of  these  results  is  pursued. 

For  this  research  Nm-lg  =  512  pixels,  and  n  is  a  power  of  two  ranging  between  1 
and  512.  Norig  must  be  divisible  by  n  so  that  Nnew  is  an  integer,  because  the  number 
of  pixels  cannot  be  fractional.  This  method  yields  ten  differently  pixelated  phase 
screens  (each  with  250  realizations)  ranging  from  lxl  pixels  to  512  x  512  pixels.  An 
example  of  one  of  the  down-sampled  screens  is  shown  in  Fig.  4.3,  and  the  original 
512  x  512  pixel  screen  is  shown  in  Fig.  4.1. 
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4-3  Finding  Phase  Variance 

The  next  step  is  to  find  the  phase  variance  across  the  250  realizations  of  each 
of  the  ten  levels  of  pixelated  phase  screens.  This  value  is  obtained  by  calculating  the 
variance  of  the  Zernike  coefficient  a%  for  each  high  order  mode  represented  across  all 
250  phase  screen  realizations.  The  number  of  modes  required  for  a  given  threshold 
of  phase  standard  deviation  desired  to  be  present  in  the  emulated  atmosphere  can 
be  found  from  Table  3.2.  Once  the  number  of  modes  is  determined,  the  Zernike 
coefficient  for  each  mode  can  be  found  for  an  individual  phase  screen  realization 
with  Eq.  (2.24).  As  shown  in  Eq.  (3.26),  the  phase  variance  of  the  ith  Zernike  mode 
is  found  by  calculating  the  statistical  variance  of  the  Zernike  coefficients  over  every 
phase  screen  realization  for  that  mode.  The  total  phase  variance  present  now  becomes 
the  sum  of  the  variances  for  each  mode.  The  Matlab®  listing  in  Appendix  C  computes 
this  phase  variance  for  each  of  the  ten  levels  of  pixelation  for  i  —  4  to  i  —  139  Zernike 
modes. 

The  simulation  analysis  was  limited  to  the  97%  threshold  level  of  total  phase 
standard  deviation  present  in  the  emulated  turbulence  due  to  the  memory  limitations 
in  Matlab®  .  In  order  for  the  listing  in  Appendix  C  to  run  in  a  reasonable  amount 
of  time,  the  Zernike  polynomials  had  to  be  pre-computed.  This  created  a  variable  in 
Matlab®  that  was  512  x  512  x  i  in  size,  where  i  is  the  number  of  modes  used.  This 
variable  became  too  large  for  Matlab®  to  keep  stored  and  continue  running  when 
the  number  of  modes  pre-computed  reached  the  98%  threshold  level  in  Table  3.2.  It 
was  possible  to  find  these  values  in  Matlab®  by  computing  each  Zernike  polynomial 
inside  of  a  loop  and  then  clearing  the  variable;  however,  the  execution  time  for  this  was 
infeasible.  It  would  require  over  two  weeks  to  run  the  program  for  all  ten  pixelation 
levels  out  to  the  99%  threshold  value  of  total  phase  standard  deviation  present.  A 
supercomputer  could  have  been  accessed  to  speed  up  this  process;  however,  the  results 
obtained  from  the  lower  threshold  levels  were  sufficient  to  forgo  pursuing  the  phase 
variance  values  for  the  98  and  99%  threshold  levels. 
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Tabic  4.1:  Simulation  results  for  the  number  of  pixels  required  to  reach  certain  threshold 

levels  of  total  phase  standard  deviation  present  after  using  a  discrete  number  of  Zernike 
modes. 


Percentage  of  total 
a<t>PTR  present 

Finite  mode 

^threshold 

j/  of  modes 
required 

Finite  pixel 

%threshold 

?7max 

#  of  pixels 
required 

99% 

99.5% 

457 

99.5% 

N/A 

N/A 

98% 

99.0% 

218 

99.0% 

N/A 

N/A 

97% 

98.5% 

139 

98.5% 

40 

25.47 

96% 

98.0% 

101 

98.0% 

33 

21.01 

95% 

97.5% 

78 

97.5% 

28 

17.83 

90% 

94.9% 

35 

94.9% 

17 

10.82 

83% 

91.1% 

20 

91.1% 

11 

7.00 

4-4  Results 

The  final  step  is  to  plot  the  results.  These  plots  are  the  result  of  the  Mat  lab® 
code  shown  in  Appendix  A.  Figure  4.4  displays  the  theoretical  and  simulation  plots 
of  normalized  phase  standard  deviation.  This  is  shown  as  a  function  of  the  number 
of  pixels  across  a  device  as  a  result  of  having  only  95%  of  the  total  phase  standard 
deviation  present.  This  is  due  to  finite  Zernike  modes  and  a  limited  spatial  frequency. 
This  phase  variance  is  normalized  because  the  D/r 0  =  5  value  was  easy  to  remove 
from  the  overall  phase  variance  after  it  was  computed. 

Figure  4.4  shows  that  the  simulation  results  match  closely  with  the  theoretical 
values  but  are  not  quite  the  same.  This  discrepancy  is  discussed  in  Chapter  V. 
There  are  only  ten  data  points  from  which  to  plot  the  simulation  results  due  to  only 
obtaining  the  phase  standard  deviation  values  at  ten  pixelation  levels  of  the  phase 
screen  realizations.  The  interp  command  in  Matlab®  is  used  to  create  the  continuous 
simulation  plot  line  in  the  figure.  This  command  conducts  a  piecewise  cubic  Hermite 
interpolation  of  the  ten  simulation  data  points.  This  interpolation  allows  for  an 
estimation  of  the  plot  values  to  be  obtained  for  every  value  of  77.  Figure  4.4  only 
shows  the  results  for  the  overall  95%  threshold  case;  however,  the  other  cases  from 
the  theoretical  analysis  shown  in  Table  3.3  are  presented  for  the  simulation  analysis 
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in  Table  4.1.  As  previously  stated,  the  98  and  99%  cases  could  not  be  evaluated 
in  simulation  due  to  the  limitations  of  the  memory  in  Matlab®  .  These  simulation 
results  are  compared  to  the  theoretical  analysis  results  in  Chapter  V. 

4-4-1  Minimum  Possible  Atmospheric  Coherence  Length  in  Simulation. 

The  example  in  Section  3. 2. 3. 4  describing  the  minimum  possible  value  of  atmo¬ 
spheric  coherence  diameter  that  a  device  can  produce  is  revisited.  This  allows  for  a 
comparison  between  the  simulation  and  theoretical  results.  In  this  example,  a  7  x  7 
pixel  device  is  theoretically  able  to  represent  96%  of  the  total  PTR  phase  standard 
deviation.  However,  this  device  can  only  represent  83%  of  the  total  PTR  phase  stan¬ 
dard  deviation  according  to  the  simulation  results  in  Table  4.1.  Therefore,  seventeen 
percent  of  the  total  is  missing  when  a  device  with  only  seven  pixels  across  is  used  in 
simulation.  The  minimum  possible  value  for  r0  that  this  device  can  represent  is  given 
in  Eq.  (3.40).  In  the  example  for  the  simulation  results,  N  =  7  and  %miSSing  =  17%. 
Equation  (3.40)  is  shown  as  a  function  of  cr<Apr™  in  Fig.  4.5  when  an  arbitrary 

x  7  '  *  i  n-missing 

value  for  the  aperture  diameter  is  set  as  D  =  10  mm.  Figure  4.5  shows  that,  in 
simulation,  a  7  x  7  pixel  device  with  an  aperture  of  D  =  10  mm  can  represent  a  min¬ 
imum  value  for  ro  =  9.856  mm  when  the  researcher  is  willing  to  allow  1/100  waves 
of  missing  PTR  phase  standard  deviation.  An  even  smaller  value  for  romin  can  be 
obtained  if  the  researcher  can  accept  a  larger  amount  of  missing  PTR  phase  standard 
deviation.  These  results  show  that  the  maximum  amount  of  turbulence  that  a  phase 
screen  can  represent  is  dependent  on  the  requirements  of  the  researcher.  The  results 
gained  from  this  example  are  discussed  in  Chapter  V  and  compared  to  the  theoretical 
results  obtained  in  Section  3. 2. 3.4. 
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(g)  GO  (i) 

Figure  4.3:  Random  individual  realization  of  a  piston/tilt  removed  phase  screen  created 

using  the  Fourier  series  method  and  then  down-sampled  to:  (a)  256  x  256  pixels  (b)  128  x  128 
pixels  (c)  64  x  64  pixels  (d)  32  x  32  pixels  (e)  16  x  16  pixels  (f)  8  x  8  pixels  (g)  4  x  4  pixels 
(h)  2x2  pixels  (i)  lxl  pixels 


56 


Number  of  Pixels,  N 


e 

<u 

C/3 

<D 

H 

H  c 

&H  O 


tn 


g  3 
Q  % 
■g  £ 

cd 

C 

•S 

55 

D 

c/3 

o3 

-C 

Oh 


Figure  4.4:  Normalized  phase  standard  deviation  present  as  a  function  of  the  number  of 
pixels  across  a  wavefront  control  device  after  only  having  95%  of  the  total  phase  standard 
deviation  present  due  to  finite  Zernike  modes  and  a  limited  spatial  frequency. 
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[waves] 


Figure  4.5:  Minimum  atmospheric  coherence  length  as  a  function  of  missing  PTR  phase 
standard  deviation  for  simulation  results  for  a  7  x  7  pixel  device  with  an  aperture  diameter 
D  =  10  mm. 


58 


V.  Conclusions  and  Future  Work 


This  chapter  revisits  the  key  results  of  this  research  and  presents  overall  con¬ 
clusions.  Further,  it  addresses  the  main  challenges  faced  during  this  research 
and  how  well  the  initial  goals  were  satisfied.  Finally,  this  chapter  presents  some  rec¬ 
ommendations  for  extensions  to  this  research  as  well  as  future  work  that  could  be 
pursued. 

5.1  Results  and  Conclusions 

5.1.1  Limitations  of  Current  Low-cost  Wavefront  Control  Devices. 

Recall,  Chapter  If  presented  a  review  of  relevant  theory  and  literature  for  this 
research.  This  review  identified  several  key  limitations  of  segmented  wavefront  control 
devices  already  researched  and  highlighted  knowledge  gaps  that  are  bridged  in  this 
thesis.  It  was  found  that  there  are  two  main  types  of  segmented  wavefront  control 
devices,  each  with  its  own  advantages  and  disadvantages.  These  two  types,  NLC  SLMs 
and  MEMS  DMs,  are  both  compact,  inexpensive,  electronically  addressable,  highly 
reproducible,  and  suited  for  AO  research  and  turbulence  emulation  in  a  laboratory 
environment. 

NLC  SLMs  have  much  larger  arrays  of  independently  controllable  elements  (at 
least  512x512  pixels)  without  any  moving  parts.  Disadvantages  of  NLC  SLMs  are  that 
they  have  slow  relaxation  times  and,  because  of  this,  they  cannot  currently  be  used 
in  real-time  AO  systems.  NLC  SLMs  are  also  wavelength  and  polarization  dependent 
devices  and  have  digitization  limitations  with  only  128  programmable  levels.  MEMS 
DMs  have  the  advantage  of  extremely  high  frame-rates  and  can  be  used  in  real-time 
AO  systems.  Segmented  MEMS  DMs  do  not  have  interactuator  coupling,  but  are 
limited  by  diffraction  effects.  Using  lenslets  to  focus  the  incident  light  in  the  center 
of  each  mirror  segment  has  been  shown  to  mitigate  these  adverse  effects. 

This  review  did  not  reveal  any  in-depth  research  on  the  effects  of  pixelation  on 
abilities  of  NLC  SLMs  and  MEMS  DMs.  This  thesis  research  bridged  that  knowledge 
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gap  through  an  investigation  of  the  limitations  of  Nyquist  bandwidth  as  a  result  of 
pixelation  for  segmented  wavefront  control  devices. 

5.1.2  Theoretically  Emulating  the  Atmosphere. 

Chapter  III  has  examined  the  theory  of  layered  turbulence  to  uncover  the  max¬ 
imum  turbulence  that  could  be  emulated  with  only  two  phase  screens.  It  was  desired 
to  find  the  minimum  possible  atmospheric  coherence  diameter  and  maximum  possible 
Rytov  number  that  could  be  achieved  with  two  phase  screens.  This  research  uncov¬ 
ered  that  both  of  these  important  atmospheric  turbulence  parameters  were  functions 
of  the  minimum  possible  atmospheric  coherence  diameter  that  a  single  segmented 
wavefront  control  device  could  produce.  Additionally,  it  was  found  that  the  Rytov 
number  is  at  a  maximum  when  both  phase  screens  in  a  system  are  at  the  beginning  of 
the  optical  propagation  path.  This  is  difficult  in  a  laboratory  due  to  the  physical  size 
of  wavefront  control  devices.  Therefore,  Fig.  3.1  presents  a  contour  plot  identifying 
the  maximum  possible  Rytov  number  a  system  could  achieve  based  on  the  position 
of  each  phase  screen. 

Another  important  result  of  the  theoretical  analysis  was  that  minimum  atmo¬ 
spheric  coherence  diameter  is  a  function  of  Nyquist  frequency  which  is  directly  related 
to  the  number  of  pixels  across  a  wavefront  control  device.  An  expression,  given  by 
Eq.  (3.22),  was  found  for  the  relationship  between  romin  and  the  Nyquist  bandlimit 
of  a  device  at  a  single  point  on  a  device.  This  presented  the  first  of  many  challenges 
in  this  research  because  without  a  detailed  correlation  matrix  of  every  point  on  a 
wavefront  control  device,  this  relationship  does  not  have  a  practical  application. 

Phase  standard  deviation  threshold  levels  and  Zernike  modes  were  pursued  to 
find  a  more  practical  relationship  of  atmospheric  coherence  diameter  to  the  number 
of  pixels  on  a  device  over  the  entire  aperture  of  the  device.  An  expression  for  the  total 
amount  of  PTR  phase  standard  deviation  present  in  a  Zernike  mode  as  a  function  of 
normalized  spatial  frequency  was  developed  and  is  shown  in  Eq.  (3.37).  This  result 
was  confirmed  to  be  accurate  through  comparison  with  a  similar  expression  given  by 
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Roggemann  and  Welsh  [31]  that  assumed  infinite  bandwidth.  Matlab®  was  then  used 
to  perform  a  numeric  analysis  to  find  the  total  percentage  of  PTR  phase  standard 
deviation  a  device  can  represent  as  a  function  of  the  number  of  pixels  across  the  device. 
Finally,  an  example  was  presented  in  Section  3. 2. 3. 4  to  show  how  this  relationship 
can  lead  to  the  minimum  possible  atmospheric  coherence  diameter  a  single  wavefront 
control  device  can  represent.  This  representation  is  a  function  of  the  number  of  waves 
of  missing  PTR  phase  standard  deviation  a  researcher  is  willing  to  accept. 

This  roundabout  numeric  method  for  finding  an  expression  for  the  minimum 
coherence  diameter  of  a  device  is  necessary  because  it  is  not  possible  to  derive  a  direct 
analytic  expression.  The  results  from  this  numeric  analysis  can  be  substituted  into 
Eqs.  (3.2)  and  (3.6)  to  obtain  the  minimum  overall  ro  and  maximum  Rytov  number 
for  a  two  phase  screen  system.  This  can  yield  the  true  limitations  in  the  strength  of 
turbulence  a  two  phase  screen  system  can  emulate  when  using  segmented  wavefront 
control  devices  as  phase  screens.  The  results  depend  on  the  acceptable  limitations  of 
the  researcher  and  the  specific  wavefront  control  device  used. 

5.1.3  Emulating  the  Atmosphere  in  Simulation. 

Chapter  IV  provides  a  confirmation  of  the  results  obtained  theoretically  through 
Matlab®  simulation.  Random  realizations  of  PTR  phase  screens  were  created  so 
that  the  phase  standard  deviation  over  those  screens  could  be  calculated.  These 
realizations  were  then  individually  down-sampled  to  several  pixelation  levels.  An 
expression  for  the  amount  of  PTR  phase  standard  deviation  present  in  a  simulated 
device  as  a  function  of  the  number  of  pixels  across  that  device  was  then  interpolated 
from  the  simulation  results  using  Matlab®  .  The  simulation  was  unfortunately  limited 
to  the  97%  threshold  level  of  total  PTR  phase  standard  deviation  because  of  the 
limitations  in  Matlab®’s  memory  capability.  This  limitation  prevented  a  comparison 
between  the  98  and  99%  threshold  level  results  from  theory  and  simulation.  Finally, 
the  same  example  pursued  in  Chapter  111  was  revisited  for  the  simulation  results. 
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Tabic  5.1:  Comparison  of  theoretical  and  simulation  results  for  the  number  of  pixels 

required  to  reach  certain  threshold  levels  of  total  phase  standard  deviation  present  after 
using  a  discrete  number  of  Zernike  modes. 


%  of  total 
a<t>PTn  present 

Finite  mode 
&  pixel  %thresh 

#  of  modes 
required 

Pixels  required 
in  theory 

Pixels  required 
in  simulation 

99% 

99.5% 

457 

15.28 

N/A 

98% 

99.0% 

218 

10.19 

N/A 

97% 

98.5% 

139 

8.28 

25.47 

96% 

98.0% 

101 

7.00 

21.01 

95% 

97.5% 

78 

6.37 

17.83 

90% 

94.9% 

35 

4.46 

10.82 

83% 

91.1% 

20 

3.82 

7.00 

5.1.4  Discrepancies  Between  Theoretical  and  Simulation  Analysis. 

The  numerical  evaluation  and  Mat  lab®  simulations  disagree  about  the  amount 
of  PTR  phase  standard  deviation  present  in  a  device.  The  plot  of  this  bandlimited 
PTR  phase  standard  deviation  shown  in  Fig.  4.4  shows  a  discrepancy  between  the  two 
different  methods.  This  discrepancy  can  also  be  visualized  in  the  last  two  columns  of 
Table  5.1  which  is  a  combination  of  the  results  from  Tables  3.3  and  4.1.  This  may 
appear  to  be  large  at  low  spatial  frequencies,  but  it  can  be  explained  quite  simply. 
Figure  4.4  shows  that  a  7  x  7  pixel  device  can  theoretically  represent  96%  of  the  total 
PTR  phase  standard  deviation  present  in  the  atmosphere.  This  is  in  contrast  to  the 
simulation,  where  it  would  take  a  21  x  21  pixel  device  to  represent  the  same  amount 
of  atmospheric  turbulence.  Table  5.1  further  shows  that,  in  simulation,  a  7  x  7  pixel 
device  is  only  capable  of  representing  83%  of  the  total  PTR  phase  standard  deviation 
present  in  the  atmosphere. 

Part  of  the  discrepancy  between  the  theoretical  and  simulation  results  can  be 
justified  by  looking  at  how  each  method  filters  out  the  threshold  levels.  The  theoretical 
analysis  numerically  evaluates  an  integral  that  goes  from  zero  to  infinity  when  all  the 
PTR  phase  standard  deviation  is  present.  This  analysis  simply  cuts  off  the  integral 
at  a  maximum  normalized  spatial  frequency  r/max  to  find  the  amount  present  up  to  a 
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certain  number  of  pixels.  This  acts  as  a  rectangle  filter  function  on  the  equation  and 
is  shown  by 

pr]  max 

^PTRpres  OC  /  f(v)d'T]  (5.1) 

Jo 

oc  [  f(rj) rect  (— - dij,  (5.2) 

Jo  V^max 

where  f(rj)  is  a  function  of  normalized  spatial  frequency  and  the  rectangle  function 
filters  the  integral  so  that  it  is  only  evaluated  from  zero  to  r]max.  This  filter  function 
causes  the  sharp  drop  seen  in  the  theoretical  results  shown  in  Fig.  4.4.  The  simulation 
method  does  not  use  a  rectangle  function  to  filter  out  the  unwanted  phase  standard 
deviation,  ft  merely  calculates  the  phase  standard  deviation  of  random  realizations 
of  phase  screens  at  different  pixelation  levels,  ft  is  shown  in  Fig.  4.4  that  as  the 
number  of  pixels  on  a  device  increases,  the  difference  between  the  simulation  and 
theoretical  results  decreases.  This  means  that  the  discrepancy  at  low  normalized 
spatial  frequencies  is  caused  by  the  pixelation  process.  In  other  words,  the  pixelation 
process  acts  as  a  filtering  function  that  removes  more  phase  standard  deviation  than 
the  rectangle  function  described  in  Eq.  (5.1).  Therefore,  it  takes  more  pixels  to  make 
up  for  the  phase  lost  from  down-sampling  the  phase  screen. 

This  result  can  be  validated  by  considering  the  down-sampling  process.  It  is 
simply  an  averaging  procedure  of  the  phase  over  large  sections  of  the  aperture.  This 
averaging  procedure  has  adverse  affects  that  reduce  the  amount  of  PTR  phase  stan¬ 
dard  deviation  calculated  across  many  phase  screen  realizations.  It  is  clear  that  an 
increase  in  the  number  of  pixels  that  are  averaged  together  leads  to  an  increase  in  the 
discrepancy  between  the  theoretical  and  simulation  results  of  PTR  phase  standard 
deviation.  The  derivation  of  the  actual  shape  of  the  filter  function  from  pixelation  of 
the  phase  screens  is  left  for  future  research.  However,  a  visual  approximation  of  these 
filter  functions  is  shown  in  Fig.  5.1.  Figure  5.1  explains  part  of  the  reason  for  why 
pixelation  results  are  different  than  the  theoretical  results.  This  process  also  explains 
some  of  the  discrepancy  between  the  examples  at  the  end  of  Chapters  III  and  IV 
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Figure  5.1:  Estimation  of  the  filter  functions  applied  in  theory  and  in  simulation  that 

result  in  discrepancies  between  the  results  of  the  two  different  methods  at  low  spatial  fre¬ 
quencies. 


Table  5.2:  Comparison  of  theoretical  and  simulation  results  for  romin  of  a  7  x  7  pixel 

wavefront  control  device. 


Analysis 

Number  of 

rn  .  for  D  =  10  mm 

'-'min 

^/romin 

method 

pixels 

and  1/100  waves  of  oy,prrR 

1  ^missing 

Theory 

7x7 

1.736  mm 

5.759 

Simulation 

7x7 

9.856  mm 

1.015 

whose  results  are  shown  in  Table  5.2.  The  example  shows  that  a  7  x  7  pixel  device 
can  theoretically  produce  a  minimum  value  of  atmospheric  coherence  length  that  is 
over  five  times  smaller  than  the  same  device  can  produce  in  simulation. 

Another  explanation  for  part  of  the  discrepancy  between  the  theoretical  results 
and  the  simulation  results  is  the  sharp  edges  of  the  pixels.  Sharp  variations  between 
values  at  the  edges  of  two  pixels  cause  fluctuations  in  frequency  space.  The  edges 
of  the  pixels  act  as  rectangle  functions  that  become  sine  functions  when  transformed 
into  the  frequency  domain.  The  cutoff  frequency  for  each  pixel  is  twice  that  of  the 
theoretical  cutoff  Fmax  given  by  Eq.  (3.3).  This  leads  to  the  effect  shown  in  Fig.  5.2 
in  which  the  sine  function  causes  some  of  the  desired  phase  variance  to  be  missed  and 
includes  some  unwanted  phase  variance  from  higher  spatial  frequencies,  ft  is  unclear 
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Figure  5.2:  Discrepancy  in  theoretical  vs.  simulation  results  from  sharp  pixel  edges  on 

a  wavefront  control  device. 


how  much  is  truly  added  or  missing  as  a  result  of  this  phenomenon.  This  effect 
could  be  included  to  make  the  theoretical  approach  more  realistic;  however,  it  would 
require  knowledge  of  the  correlation  matrix  between  all  of  the  pixels  on  a  wavefront 
control  device.  This  is  left  as  a  topic  for  future  research.  Figure  5.2  shows  that  the 
phase  variance  present  in  a  pixelated  phase  screen  is  not  entirely  representative  of  the 
atmospheric  turbulence  being  emulated. 

The  expected  results  in  an  actual  laboratory  should  be  closer  to  those  of  the 
simulation  than  the  theory  developed  in  this  research.  This  is  due  to  down-sampling 
being  required  when  a  phase  screen  is  loaded  onto  a  wavefront  control  device  so  that  it 
matches  the  number  of  pixels  on  that  device.  Therefore,  the  pixelation  filtering  process 
also  takes  place  in  a  laboratory  environment  so  that  phase  screens  can  be  loaded 
onto  wavefront  control  devices.  The  number  of  pixels  across  a  device  directly  affects 
the  amount  of  phase  that  is  lost  in  this  process.  A  512  x  512  NLC  SLM  would  not 
experience  nearly  as  much  of  a  loss  in  phase  as  a  standard  25  to  144-mirror  segmented 
MEMS  DM.  However,  the  theory  and  simulation  results  agree  for  large  N.  According 
to  Fig.  4.4  a  standard  NLC  SLM  should  be  able  to  produce  approximately  the  same 
minimum  atmospheric  coherence  length  in  theory  as  in  simulation.  Therefore,  the 
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actual  device  in  a  laboratory  should  also  meet  these  expectations.  A  MEMS  DM 
with  only  a  few  mirror  segments  is  limited  by  the  down-sampling  procedure  used  to 
map  a  phase  screen  to  the  device.  It  is  possible  that  a  more  complex  method  of  down- 
sampling  phase  screens  would  yield  simulation  results  closer  to  those  of  the  theory. 
These  filter  functions  explain  the  difference  between  the  theoretical  and  simulation 
results  in  this  research. 

5.1.5  Practicality  of  Results. 

The  results  of  this  thesis  can  be  brought  together  in  one  simple  plot  for  the 
minimum  number  of  pixels  across  a  wavefront  control  device  for  a  given  amount  of 
missing  RMS  PTR  phase  standard  deviation.  This  plot  is  representative  of  the  worst 
possible  turbulence  conditions  and  is  shown  in  Fig.  5.3.  Figure  5.3  was  created  from 
the  simulation  results  for  the  97%  case  using  139  Zernike  modes.  The  simulation 
results  were  used  because  they  are  more  realistic  than  the  theoretical  results.  This 
plot  was  limited  to  only  139  Zernike  modes  due  to  the  Mat  lab®  memory  limitations 
explained  in  Section  4.3.  Figure  5.3  was  created  using  Eq.  (3.39)  to  find  the  total 
°>ptr  .  .  for  a  variety  of  D/r0  values.  The  values  for  %miSSing  were  converted  to 
total  pixels  across  a  wavefront  control  device  through  a  plot  similar  to  that  found  in 
Fig.  4.4. 

The  practical  application  of  these  results  can  be  shown  with  an  example.  A 
given  laboratory  setup  uses  a  laser  operating  at  a  wavelength  of  A  =  500  nm  and 
has  a  total  propagation  distance  of  L  =  2  m.  The  experiment  in  question  requires 
two  wavefront  control  devices  each  with  an  aperture  diameter  of  D  =  10  mm.  The 
researcher  wants  to  know  the  minimum  number  of  pixels  a  wavefront  control  device 
must  have  to  represent  a  value  of  D/r0  =  10  and  only  be  missing  1/20  or  1/100  waves 
of  RMS  PTR  phase  standard  deviation.  According  to  Fig  5.3,  a  7  x  7  pixel  device  is 
necessary  to  represent  D/tq  =  10  to  within  1/20  waves  and  a  30  x  30  pixel  device  is 
required  to  represent  D/r0  =  10  to  within  1/100  waves  of  RMS  PTR  phase  standard 
deviation.  The  researcher  also  wants  to  know  the  maximum  Rytov  number  this  setup 
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Figure  5.3:  Total  RMS  PTR  phase  standard  deviation  missing  a versus  the 

1  ‘‘'missing 

minimum  number  of  pixels  across  a  wavefront  control  device  required  to  represent  given 
values  of  D/r q. 
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can  emulate.  The  desired  D/tq  =  10  constraint  yields  r0min  =  1  mm  for  D  —  10  mm. 
To  achieve  the  maximum  possible  Rytov  number,  both  wavefront  control  devices  are 
placed  at  the  very  beginning  of  the  optical  path.  This  yields  a  slightly  unrealistic 
value  of  z  =  0  for  the  position  of  both  phase  screens.  The  resulting  value  for  the 
maximum  possible  Rytov  number  for  this  scenario  when  all  of  these  parameters  are 
substituted  into  Eq.  (3.6)  is  cr^max  =  0.646. 

In  conclusion,  NLC  SLMs  are  rarely  adversely  affected  by  their  Nyquist  ban- 
dlimit,  but  segmented  MEMS  DMs  do  experience  limitations  in  the  amount  of  atmo¬ 
spheric  turbulence  they  can  reproduce  as  a  result  of  their  number  of  mirror  segments. 
However,  both  devices  are  capable  of  reproducing  the  total  desired  PTR  phase  stan¬ 
dard  deviation  with  high  levels  of  fidelity.  The  requirements  of  each  individual  re¬ 
search  experiment  should  determine  if  the  abilities  of  a  particular  device  are  sufficient. 

These  results  were  much  more  promising  than  initially  anticipated.  It  was  ex¬ 
pected  that  there  would  be  severe  limitations  in  the  abilities  of  all  segmented  devices. 
This  was  not  the  case,  and  it  turns  out  that  only  segmented  MEMS  DMs  with  very 
few  mirror  segments  are  affected  by  the  Nyquist  bandlimit  of  the  device.  The  other 
advantages  and  disadvantages  of  both  devices  should  be  weighed  when  considering 
which  type  of  wavefront  control  device  is  best  suited  for  a  laboratory  experiment. 

5.2  Future  Work 

This  research  is  complete,  but  there  are  many  avenues  that  may  be  investi¬ 
gated  to  further  this  investigation.  First,  the  results  of  the  theoretical  and  simulation 
analysis  presented  can  be  confirmed  experimentally.  Time  did  not  allow  for  any  exper¬ 
iments;  however,  such  an  investigation  should  be  carried  out  for  a  variety  of  devices 
with  different  Nyquist  bandlimits  to  verify  the  effects  of  pixelation  on  a  phase  screen 
device. 

Also,  an  investigation  into  the  autocorrelation  matrix  between  individual  points 
over  an  aperture  can  be  pursued.  A  better  understanding  of  how  single  points  of  an 
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aperture  are  related  will  yield  an  analytic  expression  for  the  minimum  atmospheric 
coherence  diameter  as  a  function  of  Nyquist  bandlimit  for  a  wavefront  control  device. 
An  analytic  approach  would  not  have  the  limitations  from  finite  Zernike  modes  that 
were  revealed  in  the  numerical  analysis  in  this  research. 

Another  extension  of  this  research  is  the  investigation  of  how  different  quantities 
of  phase  screens  affect  the  overall  atmospheric  properties  that  can  be  represented  in 
a  laboratory.  Along  the  same  line  of  research,  the  best  locations  for  phase  screens  in 
a  laboratory  can  be  determined  for  different  types  of  optical  sources.  This  research 
only  looked  at  a  plane  wave,  but  a  spherical  wave  would  be  interesting  as  well. 

New  designs  of  segmented  wavefront  control  devices  could  also  be  researched. 
The  pixels  on  NLC  SLMs  and  MEMS  DMs  have  sharp  edges.  These  edges  impart 
extraneous  high  frequencies  into  the  phase  screen,  ft  is  possible  that  a  spatial  filter 
inside  a  4-f  system  could  make  a  pixelated  system  appear  continuous  and  mitigate 
these  effects. 

The  simulation  analysis  can  also  be  extended  along  several  paths.  A  supercom¬ 
puter  can  be  utilized  to  carry  out  the  simulation  results  to  higher  threshold  levels. 
Also,  an  investigation  could  be  mounted  to  fully  understand  how  the  pixelation  process 
filters  the  phase  standard  deviation  results.  This  could  lead  into  work  on  improving 
the  down-sampling  procedure  to  mitigate  the  adverse  effects  down-sampling  has  on 
the  simulation  results.  A  more  complex  method  of  pixelating  a  phase  screen  could 
yield  better  simulation  results. 

Finally,  future  work  could  be  done  to  compare  the  capabilities  of  segmented 
wavefront  control  devices  with  those  of  the  more  traditional  atmospheric  phase  wheels. 
An  investigation  in  how  wavefront  sensors  are  affected  by  segmented  wavefront  control 
devices  would  also  be  of  interest.  These  future  investigations  could  truly  uncover  the 
full  range  of  limitations  and  capabilities  of  segmented  wavefront  control  devices. 
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Appendix  A.  Matlab  Code  to  Compute  Theoretical  and  Simulation 

Results 

Listing  A.l:  Matlab®  file  to  compute  theoretical  and  simulation  results. 

7.  Theoretical  and  Simulation  Analysis 
clc  ; 

format  compact; 
warning  off  all ; 

T  =  true  ; 

F  =  false  ; 

7,  Set  Flags 

include_DrO  =  F ; 


/  Define  user  changeable  variables 

eta.max  =  le4;  %  maximum  extent  of  normalized  spatial  freq 
D_over_rO  =  40;  %  value  for  D/rO  if  it  is  included 


7. 


7. 


Uncomment  the  percent  threshold  level  to 
of  phase 

standard  deviation  present  after  finite 

7,  mode_percent_pres_THRESH  =  0.911; 
of  waves  present 

7,  mode_percent_pres_THRESH  =  0.949; 
of  waves  present 

mode_percent_pres_THRESH  =  0.975;  / 

waves  present 

7,  mode_percent_pres_THRESH  =  0.980; 
of  waves  present 

7,  mode_percent_pres_THRESH  =  0.985; 
of  waves  present 

7,  mode_percent_pres_THRESH  =  0.990; 
of  waves  present 

7,  mode_percent_pres_THRESH  =  0.995; 
of  waves  present 


be  used  for  the  amount 

modes  are  considered 
7.  83 7»  total  percentage 

7.  9  0 7.  total  percentage 

95  7»  total  percentage  of 

7.  96 7o  total  percentage 

7.  97 7o  total  percentage 

7.  98 7«  total  percentage 

7.  99 7»  total  percentage 


7.  ALSO  CHANGE: 

7.  SIM_PTR_stdev_pres_after_pixel_percent(Y,:) 

7.  PLOURDE_percent_pres_simulat  ion_pos_XX 

7.  PL0URDE_per  cent  _pr e  s_  s  imulat  ion_value_XX 

7,  PL0URDE_actual_pres_s imulat  ion_value_XX 

7.  cor_eta_PLOURDE_percent_pres_simulation_XX 

7.  c  or  _Z_mode  _PL0URDE_per  cent  _pr  e  s  _mode_equat  ion_XX 

7.  in  LINES  447  -  452 

7.  to  correct  value  for  XX  and  Y 


7.  Set  the  percentage  of  waves  allowed  missing  from  finite  number 
of  pixels  and  modes 

pixel_percent_pres_THRESH  =  mode_percent_pres_THRESH ; 
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7.  Set  plot  limits  for  later 
xmin_mod.es  =  4; 
xmax.modesl  =  200; 
xmax_modes2  =  500; 

xmax.eta  =  30 ; 

xmax.pixels  =  2* xmax_eta/pi ; 

7,  Set  up  Zernike  modes  to  be  used 

7,  Load  Zernike  index  files  (containing  n,  m  ,  &  J  values  for  ... 
each  Zernike  mode 
load( ' zernike.index .mat ’ ) ; 
load( ’ zernike_index_no_trig .mat ’ ) ; 

7,  Remove  piston  and  tilt  modes  from  the  index  files 
zernike_index  (1 : 3  ,  : )  =  []  ; 

zernike_index_no_trig  (1:2,:)  =  []  ; 

7,  Set  up  Zernike  modes  to  be  used 
n  =  zernike.index  (  :  ,1)  ; 
m  =  zernike.index  (  :  ,2)  ; 

J  =  zernike.index  (  :  ,3)  ; 
nu  =  n  +  1  ; 
alpha  =  1 6 . *  nu  ; 

n_trig  =  zernike_index_no_trig  (  :  ,  1)  ; 
m_trig  =  zernike_index_no_trig  (  :  ,  2)  ; 

J_trig  =  zernike _ index_no_t rig  (:,  3)  ; 
nu_trig  =  n_trig  +  1; 
alpha_trig  =  16.*nu_trig; 

7.  Set  up  the  normalized  spatial  frequency  dummy  variable  for  ... 
integration 
eta.min  =  eps ; 

eta.numpts  =  eta.max; 

eta.space  =  (eta.max  -  et a_min ) / et a_numpt s ; 
eta  =  linspace ( eta.min ,  eta.max ,  eta.numpts) ; 

D_Fx  =  eta. /pi;  7.  spatial  frequency 

D_over_delta_pix  =  2*eta./pi;  7.  number  of  pixels  across  for 
a  given  point  in  eta 


7. 7. 7. 7. 1 1 1 1 1 1 1 7. 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 

7.  Mode  Equations  Defined  in  Roggemann  and  Welsh  % 

7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 


7.  Define  Roggemann  Mode  Equation 
for  idxl  =  1 : length ( zernike.index ) ; 

7,  Compare  to  Roggmann  phvar  &  stdev  equations 
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mode_phvar_ALL_Rogg(idxl,:)  =  4.5e-4*pi~(8/3).*(-l).~(n(.. 
idxl)  -  m(idxl)).*alpha(idxl).*... 

( ( gamma (14/3)  . *  gamma (nu (idxl)  -  11/6)). /(( gamma  (17/6) ) 
~2 . *gamma (nu ( idxl )  +  17/6))); 

90  end 

7,  Set  up  an  if  statement  to  include  D/rO  functions 
if  include_DrO 

7,  Find  the  D/rO  weighted  phase  variance  present  as  a  ... 
function  of  the  number  of  modes 

95  cum_PTR_phvar_pres_Rogg  =  D_over_r0 ~ (5/3)  . * cumsum ( .  .  . 

mode_phvar  _ALL_Rogg  )  ;  7n[rad~2] 

else 

7.  Find  the  normalized  phase  variance  present  as  a  function 
of  the  number  of  modes 

cum_PTR_phvar_pres_Rogg  =  cumsum (mode_phvar_ALL_Rogg ) . 

;  7.  [rad  ~ 2] 

end 

100 

7.  Convert  to  phase  standard  deviation  Present  in  waves 

cum_PTR_phstdev_pres_Rogg  =  sqrt ( cum_PTR_phvar_pres_Rogg ) . / (2* 
pi ) ;  i [waves ] 

7.  Find  Missing  phase  standard  deviation  in  waves 
105  cum_PTR_phstdev_miss_Rogg  =  max ( cum_PTR_phstdev_pres_Rogg)  -  . 

cum_PTR_phstdev_pres_Rogg  ;  7,  [waves] 

7.  Find  the  percentage  of  total  PTR  phase  standard  deviation  ... 
present 

cum_PTR_phstdev_pres_Rogg_per cent  =  cum_PTR_phstdev_pres_Rogg . 

. /max ( cum_PTR_phstdev_pres_Rogg ) ; 
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115  7«  Define  the  constants  to  be  used  in  the  equations  below 
C_theory  =  ( 0 . 2073/0 . 423 ) *2 ~ ( - 1 1 /3) * sqrt (pi ) ; 

7.  Define  Plourde  Mode  Equation 

for  idxl  =  1 : length ( zernike_ index ) ; 

120 

7.  Define  the  analytic  phvar  &  stdev  equations 

mode_phvar_wh_theory_ALL ( idxl  ,  : )  =  C_theory .* alpha (..  . 

idxl)  .*(( gamma (nu(idxl)  -  11/6)  .*  gamma (14/6))  ./(... 
gamma (nu(idxl)  +  17/6).*  gamma (17/6))); 

end 

125  7.  Set  up  an  if  statement  to  include  D/rO  functions 
if  include_Dr0 
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‘/,  Find  the  D/rO  weighted  phase  variance  present  as  a  ... 
function  of  the  number  of  modes 
cum_PTR_phvar_pres_theory  =  D_over_rO  “  (5/3)  .*cumsum(.. 
mode_phvar_wh_theory  _ALL  )  ;  °/n[rad~2] 

else 

°/«  F ind  the  normalized  phase  variance  present  as 
of  the  number  of  modes 

cum_PTR_phvar_pres_theory  =  cumsum(... 

mode_phvar_wh_theory_ALL  )  ;  '/,  [rad 

end 

Convert  to  phase  standard  deviation  Present  in  waves 

cum_PTR_phstdev_pres_PLOURDE  =  sqrt ( cum_PTR_phvar_pres_theory ) 
. / (2*pi) ;  % [waves] 

Convert  to  percentage  of  total  phase  standard  deviation  present 
cum_PTR_phstdev_pres_PLOURDE_percent  =  ... 
cum_PTR_phstdev_pres_PLOURDE . /max ( . .  . 
cum_PTR_phstdev_pr es_PLOURDE ) ; 

7,  Find  positions  of  and  modes  associated  with  the  various  ... 
percentage  threshold  levels 

PLOURDE_percent_pres_mode_equat ion_pos_99  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.995)); 
PL0URDE_percent_pres_mode_equation_value_99  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_99) ; 
PL0URDE_actual_pres_mode_equation_value_99  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PL0URDE_percent_pres_mode_equation_pos_99) ; 
cor_Z_mode_PL0URDE_percent_pres_mode_equation_99  =  J(... 
PL0URDE_percent_pres_mode_equation_pos_99) ; 

PLOURDE_percent_pres_mode_equat ion_pos_98  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.990)); 
PL0URDE_percent_pres_mode_equation_value_98  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_98) ; 
PL0URDE_actual_pres_mode_equation_value_98  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PL0URDE_percent_pres_mode_equation_pos_98) ; 
cor_Z_mode_PL0URDE_percent_pres_mode_equation_98  =  J(... 
PL0URDE_percent_pres_mode_equation_pos_98) ; 

PLOURDE_percent_pres_mode_equat ion_pos_97  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.985)); 
PL0URDE_percent_pres_mode_equation_value_97  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_97 ) ; 
PL0URDE_actual_pres_mode_equation_value_97  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PL0URDE_percent_pres_mode_equation_pos_97 ) ; 


a  function 

-2] 
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c  or _Z_mode  _PLOURDE_per cent _pr e  s  _mode_equat ion_  97  = 
PLOURDE_percent_pres_mode_equat ion_pos_97)  ; 

PLOURDE_percent_pres_mode_equat ion_pos_96  =  min ( f ind ( . .  . 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.980)); 
PL0URDE_percent_pres_mode_equation_value_96  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_96) ; 
PL0URDE_actual_pres_mode_equation_value_96  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PL0URDE_percent_pres_mode_equation_pos_96) ; 
c or _Z_mode_PLOURDE_per cent _pr e s _mode_equat ion_ 96  =  J(... 
PLOURDE_percent_pres_mode_equat ion_pos_96) ; 

PLOURDE_percent_pres_mode_equat ion_pos_95  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.975)); 
PL0URDE_percent_pres_mode_equation_value_95  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_95) ; 
PL0URDE_actual_pres_mode_equation_value_95  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PLOURDE_percent_pres_mode_equat ion_pos_95) ; 
c or _Z_mode _PL0URDE_per cent _pr e s _mode_equat ion_ 95  =  J(... 
PLOURDE_percent_pres_mode_equat ion_pos_95) ; 

PLOURDE_percent_pres_mode_equat ion_pos_90  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.949)); 
PL0URDE_percent_pres_mode_equation_value_90  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (.  .  . 
PL0URDE_percent_pres_mode_equation_pos_90) ; 
PL0URDE_actual_pres_mode_equation_value_90  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 

PL0URDE_per cent _pres_mode_equation_pos_90) ; 
c or _Z_mode _PL0URDE_per cent _pr e s _mode_equat ion_ 90  =  J(... 
PLOURDE_percent_pres_mode_equat ion_pos_90) ; 

PLOURDE_percent_pres_mode_equat ion_pos_83  =  min(find(... 

cum_PTR_phstdev_pres_PLOURDE_percent  >=  0.911)); 
PL0URDE_percent_pres_mode_equation_value_83  =  ... 
cum_PTR_phstdev_pres_PLOURDE_percent  (. . . 
PL0URDE_percent_pres_mode_equation_pos_83) ; 
PL0URDE_actual_pres_mode_equation_value_83  =  ... 
cum_PTR_phstdev_pres_PLOURDE  (  . .  . 
PLOURDE_percent_pres_mode_equat ion_pos_83) ; 
c or _Z_mode _PL0URDE_per cent _pr e s _mode_equat ion_83  =  J(... 
PLOURDE_percent_pres_mode_equat ion_pos_83) ; 

7.  Plot  Plourde  mode  equation  vs  Rogg  mode  equations  to  verify  they... 

are  the  same 
figure ( 1 )  ;  elf  ; 
hold  on; 
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plot(J,  cum_PTR_phstdev_pres_Rogg_percent  ,  ’k-’,  ’ linewidth ’ . . 

,  9)  ; 

plot ( J  ,  cum_PTR_phstdev_pres_PLOURDE_percent  ,  1  c-  -  ’  ,  ’  . .  . 

linewidth ’  ,  3)  ; 
hold  of  f  ; 

h_l  =  legend ( 1  Normalized  \ s igma_ {\ phi _ {PTR }}  from  Roggemann  & 
Welsh  ’  ,  .  .  . 

’Normalized  \ s igma_ {\ phi_ {PTR }}  from  thesis  ... 
derivation’,  ’location’,  ’E’); 

7,  h_t  =  title  (’\bf  Percentage  of  Normalized  Phase  Std  ... 

Deviation  Present  as  a  Function  of  Zernike  Mode  Index  \rm’) 

/  ’Total  Phase  Std  Deviation  Expressed  as  a  Fraction  . 

of  PTR  ’  ;  .  .  . 

70  ’after  including  the  first  \iti\rm  Zernike  modes  in 

representing  the  atmosphere  ’ ;  ’  ’  ;  .  .  . 

70  [’For  error  from  having  finite  modes  =  ’  num2str... 

( 1 00*  mode  _per  cent  _pres_THRESH  ,’ 7,6 . 4g  ’  )  ’7,  of  total  Phase  Std  .. 

Deviation  Present  ’]})  ; 
ax_left  =  gca; 

xl im ( [ xmin.mode s  xmax.mode s 1 ] )  ;  yl im ( [min ( . .  . 
cum_PTR_phstdev_pres_Rogg_percent )  max ( . .  . 
cum_PTR_phstdev_pres_Rogg_percent ) ] ) ; 
h_yl  =  ylabel ({ ’ Percentage  of  Total  RMS  PTR’;  ’Phase  ... 
Standard  Deviation  Present’}); 
ax_right  =  axes (’ Posit  ion ’, get ( ax_left  ,’ Posit  ion . 
YAxisLocation  ’  ,  ’ right  ’  ,  ’Color’  ,  ’none’  ,  ’YColor’  ,  ’k’)  ; 
xl im ( [ xmin.mode s  xmax.mode s 1 ] )  ;  yl im ( [min ( . .  . 
cum_PTR_phstdev_pres_Rogg)  max ( . .  . 
cum_PTR_phstdev_pres_Rogg) ] ) ; 
h_x2  =  xlabelC’Mode  Index,  \iti\rm’); 

h_y2  =  ylabel ({’ Amount  of  RMS  PTR’;  ’Phase  Standard  ... 
Deviation  Present’;  ’[waves]’}); 
set(gca,  ’FontName  ’  ,  ’Times  New  Roman’); 

set([h_x2,  h_yl,  h_y2,  h_l],  ’FontName’  ,  ’Times  New  Roman’. 

,  ’ f ontsize  ’  ,  16)  ; 

7.  Plot  values  of  percent  threshold  levels  against  the  mode  index  . 

on  the  Plourde  mode  equations  plot 
f igure (2)  ;  elf  ; 
hold  on; 

plot(J,  cum_PTR_phstdev_pres_PLOURDE_percent  ,  ’k-’  ,  ’linewidth 

’  ,  4)  ; 

plot ( J (PL0URDE_percent_pres_mode_equation_pos_83 )  ,  .  .  . 
cum_PTR_phstdev_pres_PLOURDE_percent  (.  . . 

PL0URDE_percent_pres_mode_equation_pos_83  )  ,  ’k<  ’  ,  ’  . .  . 

linewidth ’  ,  2  ,  .  .  . 

’Markersize  ’  ,  14,  ’ MarkerFaceColor  ’  ,  ’c’,  ’... 

MarkerEdgeColor ’  ,  ’k’,  ’EraseMode  ’  ,  ’none’); 

plot ( J (PL0URDE_percent_pres_mode_equation_pos_90 )  ,  .  .  . 
cum_PTR_phstdev_pres_PLOURDE_percent  (. .  . 

PL0URDE_percent_pres_mode_equation_pos_90)  ,  ’ko  ’  ,  ’  . .  . 

1 inewidth ’  ,  2  ,  .  .  . 
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’Markersize’,  14,  ’ MarkerFaceColor  ’  ,  ’  c  ’  ,  ’  .  .  . 

MarkerEdgeColor’,  ’k’,  ’EraseMode ’ ,  ’none’); 
plot(J(PL0URDE_percent_pres_mode_equation_pos_95)  ,  .  .  . 
cum_PTR_phstdev_pres_PLOURDE_percent  (  . . . 

PL0URDE_percent_pres_mode_equation_pos_95)  ,  ’  ks 1  ,  ’ . .  . 

1 inewidth ’  ,  2  ,  .  .  . 

210  ’Markersize  ’  ,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor ’ ,  ’k’,  ’EraseMode’,  ’none’); 

plot(J(PL0URDE_percent_pres_mode_equation_pos_96)  ,  .  .  . 
cum_PTR_phstdev_pres_PLOURDE_percent  (. .  . 

PL0URDE_percent_pres_mode_equation_pos_96)  ,  ’kd’  ,  ’ . .  . 

linewidth ’  ,  2  ,  .  .  . 

’Markersize ’  ,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’); 
plot(J(PL0URDE_percent_pres_mode_equation_pos_97)  ,  .  .  . 

215  cum_PTR_phstdev_pres_PLOURDE_percent  (  .  . . 

PL0URDE_percent_pres_mode_equation_pos_97 )  ,  ’  k  ~  ’  ,  ’  . .  . 

1 inewidth ’  ,  2  ,  .  .  . 

’Markersize’,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’); 
plot ( J (PL0URDE_percent_pres_mode_equation_pos_98 )  ,  .  .  . 
cum_PTR_phstdev_pres_PLOURDE_percent  (  .  . . 

PL0URDE_percent_pres_mode_equation_pos_98  )  ,  ’k> ’  ,  ’ . .  . 

linewidth ’  ,  2  ,  .  .  . 

’Markersize’,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’); 

220  plot(J(PL0URDE_percent_pres_mode_equation_pos_99)  ,  .  .  . 

cum_PTR_phstdev_pres_PLOURDE_percent  (. .  . 

PL0URDE_percent_pres_mode_equation_pos_99)  ,  ’kv’  ,  ’... 

1 inewidth ’  ,  2  ,  .  .  . 

’Markersize’,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’); 
hold  off  ; 

h_l  =  legend (’ Normalized  \ s igma_ {\ phi }  from  thesis  derivation... 

J 

9  *  *  * 

225  [’\iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_83  ,  .  .  . 
,0/o6.4g’)  ’  modes  @  ’  ... 

num2str ( 100* PLOURDE_per cent _pres_mode _equat ion_value_83  ,  .  .  . 

’ "/. 6 . 4  g  ’  )  ’7.  and  ’  ... 

230  num2str (PL0URDE_actual_pres_mode_equation_value_83  ,  .  .  . 

’  7o6 . 4g  ’  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’]  ,  .  .  . 

[ ’ \iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_90  ,  .  .  . 
’7o6.4g’)  ’  modes  @  ’  ... 

235  num2str ( 100* PL0URDE_per cent _pr es_mode_equat ion_value_90  ,  .  .  . 

’  7. 6 . 4  g  ’  )  ’  1  and  ’  ... 

num2str (PL0URDE_actual_pres_mode_equation_value_90  ,  .  .  . 

’  7o6 . 4g  ’  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’]  ,  .  .  . 

[ ’ \ it i \rm  =  ’  ... 

240  num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_95  ,  .  .  . 

’°/o6.4g’)  ’  modes  @  ’  ... 
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num2str ( 100* PLOURDE_per cent _pres _mode_equat ion_value_95  ,  .  .  . 
’7.6.4g’)  ’ 1  and  ’  ... 

num2str (PL0URDE_actual_pres_mode_equation_value_95  ,  .  .  . 

245  ’  7.6 . 4g  1  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’  ]  ,  .  .  . 

[ ’ \iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_96  ,  .  .  . 
’7o6.4g’)  ’  modes  @  ’  ... 

num2str ( 100* PL0URDE_per cent _pres_mode _equat ion_value_96  ,  .  .  . 
250  ’7.6.4 g’  )  ’7.  and  ’  ... 

num2str (PL0URDE_actual_pres_mode_equation_value_96  ,  .  .  . 

’  7o6 . 4g  ’  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’  ]  ,  .  .  . 

[ ’ \iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_97  ,  .  .  . 
255  ’7«6.4g’)  ’  modes  @  ’  ... 

num2str ( 100* PL0URDE_per cent _pr es_mode_equat ion_value_97  ,  .  .  . 

’  7.6. 4g’)  ’%  and  ’  ... 

num2str (PL0URDE_actual_pres_mode_equation_value_97  ,  .  .  . 

’  7o6 . 4g  ’  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’  ]  ,  .  .  . 

260  [’\iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_98  ,  .  .  . 
’7o6.4g’)  ’  modes  @  ’  ... 

num2str ( 100* PL0URDE_per cent _pr es_mode_equat ion_value_98  ,  .  .  . 

’  7»  6 . 4  g  ’  )  ’7.  and  ’  ... 

265  num2str (PL0URDE_actual_pres_mode_equation_value_98  ,  .  .  . 

’  7«6 . 4g  ’  )  ’  waves  of  total  \sigma_{\phi_  {PTR}}  present  ’  ]  ,  .  .  . 

[ ’ \iti\rm  =  ’  ... 

num2str ( cor_Z_mode_PL0URDE_percent_pres_mode_equation_99  ,  .  .  . 
’7o6.4g’)  ’  modes  @  ’  ... 

270  num2str ( 100* PL0URDE_per cent _pres_mode _equat ion_value_99  ,  .  .  . 

’  7«6 . 4g  ’  )  ’7.  and  ’  ... 

num2str (PL0URDE_actual_pres_mode_equation_value_99  ,  .  .  . 

’  7o6 . 4g  ’  )  ’  waves  of  total  \sigma_{\ph.i_  {PTR}}  present  ’  ]  ,  .  .  . 

’ location  ’  ,  ’ E ’ )  ; 

275  7.  h_t  =  title  (’ Percentage  of  Normalized  Phase  Std  Deviation  . 

Present  as  a  Function  of  Zernike  Mode  Index’) 

7.  ’Total  Phase  Std  Deviation  Expressed  as  a  Fraction 

of  PTR  ’  ;  .  .  . 

7.  ’after  including  the  first  \iti\rm  Zernike  modes  in 

representing  the  atmosphere ’ ;  ’  ’ ; . . . 

7.  [’For  error  from  having  finite  modes  =  ’  num2str... 

( 1 00*  mode  _per  cent  _pres_THRESH  ,’ 7«6 . 4g  ’  )  ’  7o  of  total  Phase  Std  . 

Deviation  Present  ’]})  ; 
ax_left  =  gca; 

280  xl im ( [ xmin.mode s  xmax_modes2] )  ;  yl im ( [min ( . .  . 

cum_PTR_phstdev_pres_PLOURDE_percent )  max ( . .  . 
cum_PTR_phstdev_pres_PLOURDE_percent ) ] ) ; 
h_yl  =  ylabel ({ ’ Percentage  of  Total  RMS  PTR’;  ’Phase  ... 
Standard  Deviation  Present’}); 
ax_right  =  axes (’ Pos it  ion ’, get ( ax_left  ,’ Posit  ion ’),’.. . 
YAxisLocation ’  ,  ’ right  ’  ,  ’Color’  ,  ’none’  ,  ’YColor’  ,  ’ k ’ )  ; 
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xl im ( [ xmin.mode s  xmax_modes2] )  ;  yl im ( [min ( .  .  . 
cum_PTR_phstdev_pres_PLOURDE)  max ( . .  . 
cum_PTR_phstdev_pres_PLOURDE) ] ) ; 
h_x2  =  xlabelC’Mode  Index,  \iti\rm’); 

285  h_y2  =  ylabel  ({  ’  Amount  of  RMS  PTR 1  ;  ’  Phase  Standard  ... 

Deviation  Present’;  ’[waves]’}); 
set(gca,  ’FontHame ’  ,  ’Times  New  Roman’); 

set([h_x2,  h_yl,  h_y2,  h_l],  ’FontName’  ,  ’Times  New  Roman’. 
,  ’ f ontsize  ’  ,  16)  ; 
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7.  PLOURDE  Mode  Simulation  7. 
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7.  Define  the  percentage  threshold  modes  of  interest 

thresh_modes  =  [cor_Z_mode_PL0URDE_percent_pres_mode_equation_90. . 


cor_Z_mode_PL0URDE_percent_pres_mode_equation_95 . . 

cor_Z_mode_PL0URDE_percent_pres_mode_equation_96 . . 

300  cor_Z_mode_PL0URDE_percent_pres_mode_equation_97 . . 

cor_Z_mode_PL0URDE_percent_pres_mode_equation_98 . . 

cor_Z_mode_PL0URDE_percent_pres_mode_equation_99] 

7.  Define  the  constants  to  be  used  in  the  equations  below 
305  C_sim  =  ( 0 . 2073/0 . 423 ) *2* pi *2 ~ ( -  1 1/3 )  ; 

7.  Iterate  over  each  threshold  percentage  level  of  interest 
for  idx_thresh  =  1 : length ( thre sh_mode s ) 

310  7.  Find  the  number  of  corresponding  modes  needed  for  the  ... 

zernike_index_no_trig 

N_modes_trig ( idx_thresh)  =  min ( f ind ( J_tr ig  >=  thresh.modes 
( idx_thresh ) ) ) ; 

7,  Define  Plourde  Mode  Simulation 

for  idx2  =  1 : N_modes_trig ( idx_thresh) ; 

315 

7.  Define  the  inside  of  the  numeric  integral 
if  m_trig(idx2)  ==  0; 

SIM_mode_phvar _wh_ALL ( idx2 , : )  =  C_sim.*... 

alpha_tr ig ( idx2 )  .*(besselj  (nu_trig ( idx2 )  ,  . 

eta)) .~2.*eta.“(-14/3) ; 

else 
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SIM_mode_phvar _wh_ ALL ( idx2 , : )  =  2*C_sim.*... 
alpha_t r ig ( idx2 ) .*(besselj (nu_trig ( idx2 ) , 
eta)) ,~2.*eta.~(-14/3) ; 

end 

7,  Integrate  the  above  equation  over  eta  and  find  ... 
numeric  phvar  &  stdev 

SIM_mode_cumsum_phvar_wh_ALL ( idx2 , : )  =  cumsum(... 

SIM_mode_phvar_wh_ALL (idx2,:)).*eta_space;  ... 
7.  [rad  -2] 

end 


7.  Set  up  an  if  statement  to  include  D/rO  functions 
if  include_DrO 

7.  Find  D/rO  weighted  Simulated  PTR  phvar  present  in  ... 
terms  of  spatial  frequency  after  modes  have  been  ... 
removed 

330  SIM_PTR_phvar_pres_af ter_pixel ( idx_thresh  ,  : )  =  ... 

D_over_rO  “  (5/3)  . *sum(... 

SIM_mode_cumsum_phvar_wh_ALL )  ;  7«  [rad~2] 

else 

7.  Find  normalized  Simulated  PTR  phvar  present  in  terms, 
of  spatial  frequency  after  modes  have  been  removed 
SIM_PTR_phvar_pres_af ter_pixel ( idx_thresh , : )  =  sum. 

( SIM_mode_cumsum_phvar_wh_ALL ) ;  i  [rad~2] 

end 

335 

7,  Convert  to  Simulated  PTR  Phase  Std  Dev  present  in  terms  of  .. 
spatial  frequency  after  modes  have  been  removed 
SIM_PTR_stdev_pres_af ter_pixel ( idx_thresh , : )  =  sqrt(... 

SIM_PTR_phvar_pres_after_pixel(idx_thresh  ,  :))  . / (2*pi )  ;  . 

7.  [waves] 

7,  Find  maximum  value  of  Simulated  PTR  phase  STD  Dev  present  ... 
after  pixel  have  been  removed 

340  SIM_PTR_stdev_pres_af ter_pixel_value ( idx_thresh , : )  =  max ( . . 

SIM_PTR_stdev_pres_af  ter_pixel  (  idx_thresh  ,  :  )  )  ;  7.  [waves. 

] 

7,  Convert  to  percentage  of  total  phase  standard  deviation  ... 
present  after  pixelation  effects 
SIM_PTR_stdev_pres_af ter_pixel_percent ( idx_thresh  ,  : )  =  ... 

SIM_PTR_stdev_pres_after_pixel (idx_thresh  ,  :)  . /max ( . .  . 
SIM_PTR_stdev_pres_after_pixel (idx_thresh , :) ) ; 

end 

345 

7.  Find  positions  of  and  number  of  pixels  associated  with  the  ... 
various  percentage  threshold  levels 
PL0URDE_percent_pres_simulation_pos_90  =  min(find(... 

SIM_PTR_stdev_pres_af ter_pixel_percent ( 1 , : )  >=  0.949)); 
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PL0URDE_percent_pres_simulatioii_value_90  =  ... 
SIM_PTR_stdev_pres_af ter_pixel_percent  (1  ,  ... 

PL0URDE_percent_pres_simulation_pos_90) ; 
PL0URDE_actual_pres_simulation_value_90  =  ... 
SIM_PTR_stdev_pres_af ter_pixel  (1  ,  ... 

PL0URDE_percent_pres_simulation_pos_90) ; 
cor_eta_PL0URDE_percent_pres_simulation_90  =  eta(... 

PL0URDE_percent_pres_simulation_pos_90) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_90  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_90)  ; 

PL0URDE_percent_pres_simulation_pos_95  =  min(find(... 

SIM_PTR_stdev_pres_af ter_pixel_percent (2 , : )  >=  0.975)); 

PL0URDE_percent_pres_simulation_value_95  =  ... 
SIM_PTR_stdev_pres_after_pixel_percent  (2  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_95) ; 
PL0URDE_actual_pres_simulation_value_95  =  ... 
SIM_PTR_stdev_pres_after_pixel  (2  ,  . . . 

PL0URDE_percent_pres_simulation_pos_95) ; 
cor_eta_PL0URDE_percent_pres_simulation_95  =  eta(... 

PL0URDE_percent_pres_simulation_pos_95) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_95  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_95 ) ; 

PL0URDE_percent_pres_simulation_pos_96  =  min ( f ind (  .  .  . 

SIM_PTR_stdev_pres_af ter_pixel_percent (3 , : )  >=  0.980)); 

PL0URDE_percent_pres_simulation_value_96  =  ... 
SIM_PTR_stdev_pres_after_pixel_percent  (3  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_96) ; 
PL0URDE_actual_pres_simulation_value_96  =  ... 
SIM_PTR_stdev_pres_after_pixel  (3  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_96) ; 
cor_eta_PL0URDE_percent_pres_simulation_96  =  eta(... 

PL0URDE_percent_pres_simulation_pos_96) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_96  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_96 ) ; 

PL0URDE_percent_pres_simulation_pos_97  =  min(find(... 

SIM_PTR_stdev_pres_af ter_pixel_percent (4 , : )  >=  0.985)); 

PL0URDE_percent_pres_simulation_value_97  =  ... 
SIM_PTR_stdev_pres_after_pixel_percent (4  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_97 ) ; 
PL0URDE_actual_pres_simulation_value_97  =  ... 
SIM_PTR_stdev_pres_after_pixel  (4  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_97 ) ; 
cor_eta_PL0URDE_percent_pres_simulation_97  =  eta(... 

PL0URDE_percent_pres_simulation_pos_97 ) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_97  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_97 ) ; 

PL0URDE_percent_pres_simulation_pos_98  =  min(find(... 

SIM_PTR_stdev_pres_af ter_pixel_percent (5 , : )  >=  0.990)); 
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PL0URDE_percent_pres_simulation_value_98  =  ... 
SIM_PTR_stdev_pres_after_pixel_percent  (5  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_98)  ; 
PL0URDE_actual_pres_simulation_value_98  =  ... 
SIM_PTR_stdev_pres_after_pixel  (5  ,  . . . 

PL0URDE_percent_pres_simulation_pos_98) ; 
cor_eta_PL0URDE_percent_pres_simulation_98  =  eta(... 

PL0URDE_percent_pres_simulation_pos_98) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_98  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_98) ; 

PL0URDE_percent_pres_simulation_pos_99  =  min(find(... 

SIM_PTR_stdev_pres_af ter_pixel_percent (6 , : )  >=  0.995)); 

PL0URDE_percent_pres_simulation_value_99  =  ... 
SIM_PTR_stdev_pres_after_pixel_percent  (6  ,  . . . 

PL0URDE_percent_pres_simulation_pos_99) ; 
PL0URDE_actual_pres_simulation_value_99  =  ... 
SIM_PTR_stdev_pres_after_pixel  (6  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_99) ; 
cor_eta_PL0URDE_percent_pres_simulation_99  =  eta(... 

PL0URDE_percent_pres_simulation_pos_99) ; 
cor_pixels_PLOURDE_percent_pres_simulat ion_99  =  ... 

D_over_delta_pix (PL0URDE_percent_pres_simulation_pos_99) ; 

70  Plot  results  of  theoretical  numerical  analysis  of  phase  standard... 

deviation  a  a  function  of  the  number  of  pixels  on  a  device 
f igure (3)  ;  elf  ; 
hold  on; 

plot(eta,  SIM_PTR_stdev_pres_af ter_pixel_percent (1 , : ) ,  ’k-’,  ’... 

linewidth  ’  ,  4)  ; 

plot(eta,  SIM_PTR_stdev_pres_af ter_pixel_percent (2 , : ) ,  ’k--’,  ... 

’ linewidth  ’  ,  4)  ; 

plot (eta,  SIM_PTR_stdev_pres_after_pixel_percent(4,:),  ’k- .  ’  ,  ... 

’  linewidth  ’  ,  4)  ; 

plot (eta,  SIM_PTR_stdev_pres_after_pixel_percent(6,:),  ’  k  :  ’  ,  ’  .  .  . 

linewidth ’  ,  4)  ; 

plot ( et a ( PL0URDE_per c ent _pr es _s imulation_pos_90)  ,  .  .  . 
SIM_PTR_stdev_pres_after_pixel_percent  (1  ,  .  . . 

PL0URDE_percent_pres_simulation_pos_90)  ,  ’  ko ’  ,  ’ . .  . 

1 inewidth ’  ,  2  ,  .  .  . 

’Markersize ’  ,  14,  ’ MarkerFaceColor  ’  ,  ’  c  ’  ,  ’... 

MarkerEdgeColor ’  ,  ’k1,  ’EraseMode ’  ,  ’none’); 

plot ( et a ( PL0URDE_per cent _pr e s_ s imulat ion_po s _ 95 )  ,  .  .  . 
SIM_PTR_stdev_pres_after_pixel_percent  (2  ,  .  . . 

PL0URDE_percent_pres_simulation_pos_95)  ,  ;ks ’  ,  1  . .  . 

linewidth ’  ,  2  ,  .  .  . 

’Markersize’,  14,  ’MarkerFaceColor’,  ’c’,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode ’ ,  ’none’); 
plot (eta( PL0URDE_per cent _pre s_ s imulat ion_pos_97)  ,  .  .  . 
SIM_PTR_stdev_pres_after_pixel_percent  (4  ,  . .  . 

PL0URDE_percent_pres_simulation_pos_97  )  ,  ’ k~ ’  ,  ’  . .  . 

1 inewidth ’  ,  2  ,  .  .  . 
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’Markersize’,  14,  ’ MarkerFaceColor  ’  ,  ’  c  ’  ,  ’  .  .  . 

MarkerEdgeColor ’ ,  ’k’,  ’EraseMode ’ ,  ’none’); 

plot (eta (PL0URDE_percent_pres_simulation_pos_99  )  ,  .  .  . 
SIM_PTR_stdev_pres_after_pixel_percent  (6  ,  . . . 

PL0URDE_percent_pres_simulation_pos_99)  ,  ’  kv  ’  ,  1  . .  . 

1 inewidth ’  ,  2  ,  .  .  . 

’Markersize ’  ,  14,  ’MarkerFaceColor’,  ’  c  ’  ,  ’... 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’); 
hold  off  ; 

h_l  =  legend (’ Normalized  \ s igma_ {\ phi _ {pixel s }}  as  a  function 
of  N  pixels  for  907.  of  the  total  \ s igma_ {\ phi_ {PTR}}  ’  ,  .  .  . 

’Normalized  \ s igma_ {\ phi_ {pixels }}  as  a  function  . 
of  N  pixels  for  957,  of  the  total  \sigma_{\phi_{ 
PTR}}  ’  ,  .  .  . 

’Normalized  \ s igma_ {\ phi_ {pixels }}  as  a  function  . 
of  N  pixels  for  97 7»  of  the  total  \sigma_{\phi_{ 
PTR}}  ’  ,  .  .  . 

’Normalized  \ s igma_ {\ phi_ { pixel s }}  as  a  function  . 
of  N  pixels  for  997,  of  the  total  \sigma_{\phi_{ 
PTR}}  ’  ,  .  .  . 

[’\eta  =  ’  num2str(... 

cor_eta_PL0URDE_percent_pres_simulation_90  ,  ’7.6.4g’ 

)  .  .  . 

’  &  N  =  ’  num2str  ( . .  . 

cor_pixels_PL0URDE_percent_pres_simulation_90  ,  ’... 

7.6. 4g’)  .  .  . 

’  pixels  @  ’  num2str ( 100* . .  . 

PL0URDE_percent_pres_simulation_value_90  ,  ’  7.6 . 4g  ’  )  . 

’7.  and  ’  num2str(... 

PL0URDE_actual_pres_simulation_value_90  ,  ’7.6, 4g’)  ’ 

waves  of  total  \sigma_{\phi_{PTR}}  present’],... 
[’\eta  =  ’  num2str(... 

cor_eta_PL0URDE_percent_pres_simulation_95 ,  ’.. 

7. 6 . 4  g  ’  )  .  .  . 

’  &  N  =  ’  num2str ( . .  . 

cor_pixels_PLOURDE_percent_pres_simulat ion_9  5  , 

’  7. 6 . 4  g  ’  )  .  .  . 

’  pixels  @  ’  num2str  ( 100* . .  . 

PL0URDE_per  cent  _pr  es_s  imul  at  ion_  value_95  ,  ’7.6.4 

g’)  •  •  • 

’7.  and  ’  num2str(... 

PL0URDE_actual_pres_simulation_value_95  ,  ’7.6. 4g 

’)  ’  waves  of  total  \ s igma_ {\ phi_ {PTR }}  present 

’],... 

[’\eta  =  ’  num2str(... 

cor_eta_PL0URDE_percent_pres_simulation_97  ,  ’7.6. 4g’) 

’  &  N  =  ’  num2str  (  . .  . 

cor_pixels_PL0URDE_percent_pres_simulation_97  ,  ’... 

7. 6 . 4  g  ’  )  .  .  . 
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’  pixels  @  1  num2str ( 100* . .  . 

PL0URDE_percent_pres_simulation_value_97  ,  1  7.6 . 4g  ’  )  . 


420 


7. 


425  7, 
7. 
7. 
7. 


430 


435 


’7.  and  ’  num2str(... 

PL0URDE_actual_pres_simulation_value_97  ,  5  7o 6 . 4 g  1  )  1 

waves  of  total  \ s igma_ {\ phi_ {PTR }}  present ’]  ,  .  .  . 
[’\eta  =  ’  num2str(... 

cor_eta_PL0URDE_percent_pres_simulation_99  ,  ’... 

7.6. 4g’)  .  .  . 

’  &  N  =  ’  num2str ( . .  . 

cor_pixels_PLOURDE_percent_pres_simulat ion_9  9  , 

’  7. 6 . 4  g  ’  )  .  .  . 

’  pixels  @  1  num2str ( 100* . .  . 

PL0URDE_per  cent  _pr  es  _s  imulat  ion_  value_99  ,  ’7.6.4 

g’)  •  •  • 

’7.  and  ’  num2str(... 

PL0URDE_actual_pres_simulation_value_99  ,  ’7.6. 4g 

’)  ’  waves  of  total  \  s  igma_  {\ phi  _  {PTR }}  present 

’],... 

’location’,  ’ E ’ ) ; 

h_t  =  title ({’ Percentage  of  Finite  Mode  Limited  Normalized  . 
Phase  Std  Deviation  Present ’;’ as  a  Function  of  Number  of  Pixels 
across  a  Wavefront  Control  Device  ’  ;  .  .  . 

J  J  .  J  J  .  J  J  \  . 

>  i  J  J  y 

axpos  =  get ( gca ,  ’  pos  ’  )  ; 
extent  =  get (h_t, ’extent’); 

set (gca ,’ pos ’,  [axpos (1)  axpos (2)  axpos  (3)  axpos (4)  - (4/5) 
♦extent (4) ] ) 
ax_left  =  gca; 

xlim([0  xmax.eta])  ;  ylim ( [min (min  (  . . . 

SIM_PTR_stdev_pres_after_pixel_percent ) )  max (max ( . .  . 
SIM_PTR_stdev_pres_after_pixel_percent ) ) ] ) ; 
h_xl  =  xlabel (’ Normalized  Spatial  Frequency,  \eta’); 
h_yl  =  ylabel ({’ Percentage  of  Total  RMS  PTR’;  ’Phase  ... 
Standard  Deviation  Present’}); 
ax_right  =  axes (’ Pos it  ion ’, get ( ax_left  ,’ Posit  ion ’),’.. . 

XAxisLocat ion  ’  ,  ’top’,  ’ YAxisLocation ’  ,  ’right’,  ’Color’,  ’ 
none’,  ’ YColor ’  ,  ’k’); 

xlim([0  xmax.pixels ] )  ;  ylim ( [min (min (..  . 

SIM_PTR_stdev_pres_after_pixel ) )  max (max ( . .  . 
SIM_PTR_stdev_pres_after_pixel ) ) ] ) ; 
h_x2  =  xlabel ({’ Number  of  Pixels,  N’;}); 

h_y2  =  ylabel ({’ Amount  of  RMS  PTR’;  ’Phase  Standard  ... 
Deviation  Present’;  ’[waves]’}); 
set(gca,  ’FontName’  ,  ’Times  New  Roman’); 

set([h_xl,  h_x2,  h_yl  ,  h_y2,  h_l]  ,  ’FontName’  ,  ’Times  New  . 
Roman’,  ’ fontsize ’  ,  16); 
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445  7.  MANUALLY  SET  the  FOLLOWING  VALUES  to  correspond  with  desired  ... 
threshold  level 

7,  This  is  done  for  compar ison  with  theoretical  results 
SIM_PTR_stdev_pres_af ter_pixel_per cent  =  ... 

SIM_PTR_stdev_pres_af ter_pixel_percent  (2,  :)  ; 
PLOURDE_percent_pres_simulation_pos  =  ... 

PL0URDE_percent_pres_simulation_pos_95; 
PLOURDE_percent_pres_simulation_value  =  ... 
PL0URDE_percent_pres_simulation_value_95; 

450  PLOURDE_actual_pres_simulation_value  =  ... 

PL0URDE_actual_pres_simulation_value_95; 
cor_eta_PLOURDE_percent_pres_simulation  =  ... 

cor _et a_PLOURDE_per cent  _pr es  _s imulat ion_  9  5 ; 
cor_Z_mode_PLOURDE_percent_pres_mode_equation  =  ... 
cor_Z_mode_PL0URDE_percent_pres_mode_equation_95; 

°l  Calculate  values  to  use  in  plots 

455  cor_pixels_PLOURDE_percent_pres_simulat ion  =  D_over_delta_pix ( 

PLOURDE_percent_pres_simulation_pos) ; 

°i  Load  the  necessary  files  created  from  the  computer  s  imulat  ion 
load ( ’ Thes i s\Matlab\mat  Files\HO_computer_SIM_results . mat ’ ) 

460  /  Find  the  correct  standard  deviation  values  that  match  the  number 
of  modes  used  for  the  threshold  level  of  interest 
Comput er _S IM_phs t de v_pr e s _ af t er .mode s_and_p ixel (2 : 1 1 , 2)  =  ... 

HO_norm_cum_phstdev_a_mode_ALL  (  :  ,  .  . . 
cor_Z_mode_PLOURDE_percent_pres_mode_equation -3) ; 

7,  Add  a  zero  pixel  at  ion  value  that  to  provide  more  data  points 
for  the  interpolation 

Computer_SIM_phstdev_pres_af ter_modes_and_pixel (1  ,  2)  =  0; 

465 

7,  Set  up  an  if  statement  to  include  D/rO  functions 
if  include_DrO 

7.  Find  D/rO  weighted  Comp  SIM  PTR  phstdev  present  in  terms 
of  spatial  frequency  after  modes  have  been  removed 
comp_SIM_value_pres  =  D_over_r0 ~  (5/6)  .*... 

Computer_SIM_phstdev_pres_after_modes_and_pixel . .  . 

(  :  ,  2)  ;  7.  [waves  ] 

470  else 

7.  Find  normalized  Comp  SIM  PTR  phstdev  present  in  terms  of 
spatial  frequency  after  modes  have  been  removed 
comp_SIM_value_pres  =  ... 

Computer_SIM_phstdev_pres_after_modes_and_pixel . .  . 

(  :  ,  2)  ;  7.  [waves  ] 

end 

475  7.  Create  array  of  values  so  interpolation  goes  out  to  the  correct 
number  of  data  points 
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comp_SIM_pixels  =  ... 

Computer _SIM_phstdev_pres_af ter_modes_and_pixel  (  :  ,  1)  ; 
comp_SIM_eta  =  ( pi /2) . * comp_S IM_pixel s ; 

•/.  Interpolate  the  amount  present  to  work  for  all  values  of  eta 
480  comp_SIM_value_pres_INTERP  =  int erp 1 ( comp_S IM_et a  (  1  :  1 1 )  ,  ... 

comp_SIM_value_pres(l:ll),  eta,  ’cubic’);  7o  [waves] 

7,  Convert  to  percentage  of  total  phase  standard  deviation  present 
after  pixelation  effects 

comp_SIM_percent_pres  =  comp_SIM_value_pres . /max (  .  .  . 
comp_SIM_value_pres ) ; 

comp_SIM_percent_pres_INTERP  =  comp_SIM_value_pres_INTERP . /max 
( comp_SIM_value_pres_INTERP ) ; 

485 

7.  Find  positions  of  and  number  of  pixels  associated  with  the  ... 
various  percentage  threshold  levels 
COMP_percent_pres_simulation_pos  =  min(find(... 

comp_SIM_percent_pres_INTERP  >=  pixel_percent_pres_THRESH) ) 

i 

COMP_percent_pres_simulation_value  =  ... 
comp_SIM_percent_pres_INTERP  (  . . . 
COMP_percent_pres_simulation_pos) ; 
COMP_actual_pres_simulation_value  =  comp_SIM_value_pres_INTERP 
(COMP_percent_pres_s imulat ion.pos) ; 

490  cor_eta_COMP_percent_pres_simulation  =  eta(... 

COMP_percent_pres_simulation_pos) ; 
cor_pixels_COMP_percent_pres_simulation  =  D_over_delta_pix  (  . .  . 
COMP_percent_pres_simulation_pos) ; 

7.  Plot  computer  simulation  vs.  theoretical  results  for  phase  std  . 
dev  present  after  pixelation  effects  as  a  function  of  the  ... 
number  of  pixels  on  a  device 
f igure (4)  ;  elf  ; 

495  hold  on; 

plot(eta,  SIM_PTR_stdev_pres_af ter_pixel_percent  ,  ’k-’  ,  ’. 

linewidth ’  ,  4)  ; 

plot(eta,  comp_SIM_percent_pres_INTERP  ,  ’k--’,  ’.. 

linewidth ’  ,  4)  ; 

plot ( comp_SIM_eta ,  comp_SIM_percent_pres  ,  ’ ko ’  ,  .  .  . 

’Markersize ’  ,  14,  ’MarkerFaceColor  ’  ,  ’k’,  ’.. 

MarkerEdgeColor ’ ,  ’k’,  ’EraseMode ’ ,  ’none’ 

) ; 

500  plot (eta(PLOURDE_percent_pres_s imulat ion.pos)  ,  .  .  . 

SIM_PTR_stdev_pres_after_pixel_percent  (. . . 

PLOURDE_percent_pres_simulation_pos) , ’ks’ . 
,  ’ linewidth ’  ,  2  ,  .  .  . 

’Markersize’,  14,  ’MarkerFaceColor’,  ’ c  ’  ,  ’.. 

MarkerEdgeColor’,  ’k’,  ’EraseMode’,  ’none’ 

) ; 

plot (eta(COMP_percent_pres_s imulat ion.pos)  ,  .  .  . 
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comp_SIM_percent_pres_INTERP  (. . . 

COMP_percent_pres_simulation_pos)  ,  ’  kd  ’  ,  ’  . .  . 

linewidth ’  ,  2  ,  .  .  . 

’ Markersize '  ,  14,  ’MarkerFaceColor ’  ,  ’  c  ’  ,  ’  .  .  . 

MarkerEdgeColor 1 ,  ’k’,  ’EraseMode ’ ,  ’none’... 

) ; 

hold  of  f  ; 

h_l  =  legend (’ Theoretical  normalized  \ s igma_ {\ phi }  present... 
} 

>  •  •  • 

’Cubic  interpolation  of  computer  simulated  ... 

normalized  \ s igma_ {\phi }  present  ’  ,  .  .  . 

’Actual  data  values  of  computer  simulated  ... 
normalized  \ s igma_ {\phi }  present ’  ,  .  .  . 

[’\eta  =  ’  num2str(... 

cor _eta_PLOURDE_per cent _pres_simulat ion  ,  ’  .  . . 

"/o  6 . 4  g  ’  )  .  .  . 

’  &  N  =  ’  num2str ( . .  . 

cor_pixels_PLOURDE_percent_pres_simulation  .  . . 

,  ’  °/o  6 . 4  g  ’  )  .  .  . 

’  pixels  0  ’  num2str  ( 100* . .  . 

PLOURDE_percent_pres_simulation_value  ,  ’... 

%  6 . 4  g  ’  )  .  .  . 

’  7«  and  ’  num2str(... 

PLOURDE_actual_pres_simulation_value  ,  ’7,6.4... 

g  ’  )  ’  total  phstd  present ’]  ,  .  .  . 

[’\eta  =  ’  num2str(... 

cor_eta_COMP_percent_pres_simulation  ,  ’  7,6 . 4g  ’  )  .  .  . 

’  &  N  =  ’  num2str ( . .  . 

cor_pixels_COMP_percent_pres_s  imulat  ion  ,  ’  7,6 . 4g  .  .  . 

’  )  .  .  . 

’  pixels  @  ’  num2str  ( 100* . .  . 

COMP_percent_pres_simulation_value  ,  ’7,6. 4g’)  .  .  . 

’7,  and  ’  num2str  (  C0MP_actual  _pre  s_s  imulat  ion_  value  ..  . 
,  ’  7,6 . 4g  ’  )  ’  total  phstd  present’],  ’location’... 

,  ’  E  ’  )  ; 

7,  h_t  =  title ({’ Normalized  phase  standard  deviation  ... 

present  as  a  function  of  the  number  of  pixels  across  a  ... 
wavefront  control  device  ’ ;  ’  ’  ;  .  .  . 

7.  ’after  only  having  957,  of  the  total  \... 

phi_{PTR}  present  due  to  finite  Zernike  modes  and  limited  ... 
spatial  frequency  ’  :  .  .  . 

7. 

7,  axpos  =  get  (gca  ,  ’  pos  ’  )  ; 

7,  extent  =  get  (h_t  ,  ’  extent  ’  )  ; 

7,  set(gca,’pos’,[axpos(l)  axpos  (2)  axpos  (3) 

7,  axpos  (4)  -  (5/6)  *  ext  ent  (4)  ]  ) 

ax_left  =  gca; 

xlim([0  xmax.eta] )  ;  ylim ( [min (min (  .  .  . 

SIM_PTR_stdev_pres_after_pixel_percent ) )  max (max ( . .  . 
SIM_PTR_stdev_pres_after_pixel_percent ) ) ] ) ; 
h_xl  =  xlabel (’ Normalized  Spatial  Frequency,  \eta’); 
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h_yl  =  ylabel ({’ Percentage  of  Total  RMS’;  ’Phase  ... 
Standard  Deviation  Present’}); 
ax_right  =  axes (’ Posit  ion ’, get (  ax_left  ,’ Posit  ion ’),’.. . 

XAxisLocation  ’  ,  ’top’,  ’  YAxisLocation ’  ,  ’right’,  ’Color. 
’,  ’none’,  ’ YColor  ’  ,  ’k’); 

xlim([0  xmax.pixels] )  ;  ylim ( [min (min ( . .  . 

SIM_PTR_stdev_pres_after_pixel) )  max (max (. .  . 
SIM_PTR_stdev_pres_after_pixel) )] ) ; 
h_x2  =  xlabel (’ Number  of  Pixels,  N’); 

h_y2  =  ylabel ({’ Amount  of  RMS  PTR ’ ;  ’Phase  Standard  ... 
Deviation  Present  ’  ;  ’  [waves]  ’})  ; 

set (gca ,  ’FontName’  ,  ’Times  New  Roman’); 

set([h_xl,  h_x2,  h_yl,  h_y2,  h_l],  ’FontName’  ,  ’Times  .. 
New  Roman’,  ’ f ontsize  ’  ,  16); 
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Appendix  B.  Matlab  Code  to  Create  Phase  Screens  with  the  Fourier 

Series  Method 


Listing  B.l:  Matlab®  file  to  create  phase  screens  using  the  Fourier  series  method. 

7«  Fourier  Series  Phase  Screen  Creator 
clear  all ;  clc  ; 


7.  load  path  for  locations  of  f  unct  ions  called 
addpathC’ThesisXMatlab’) ; 


7.  Define  important  parameters 

numscreens  =  250;  %  number  of  random  phase  screens  to  create 

N  =  512;  1  number  of  screen  grid  points 

D  =  2;  7,  aperture  diameter  [m] 

D_over_r0  =  5  70  strength  of  turbulence  to  be  created 

idxStr_DrO  =  sprintf  (  ’  7c02d  ’  ,  D_over_r0); 


7.  Set  the  limits  of  integration 
eta.min  =  eps ; 
eta.max  =  1 e4 ; 
eta.numpts  =le4; 

eta.space  =  (eta.max  -  et a_min ) / et a_numpt s ; 
eta  =  linspace ( eta.min ,  eta.max ,  eta_numpts); 

7.  Calculate  other  important  parameters 

rO  =  D  . /D_over_rO  ;  70  atmospheric  coherence  diameter  [m] 

deltax  =  D/N ;  70  grid  spacing 

Fx  =  eta./(D*pi);  7»  linear  spatial  frequency 

7«  Create  a  circular  aperture  of  unit  radius 
x  =  linspace (-D/2 , D/2  ,  N); 

[x  y]  =  meshgrid (x , x)  ; 

[theta  rho]  =  cart2pol (x , y ) ; 
circ_ap  =  rho  <  D/2; 
r  =  (2 . *  rho )  . /D ; 


7c  Define  the  piston  and  tilt  Zernike  mo36des  to  be  removed 
Z_piston  =  1 ; 

Z_xtilt  =  2 . *r . *  cos (theta)  ; 

Z_ytilt  =  2 . * r . * s in ( thet a ) ; 

7c  PSD  setup 

T.lo  =  0;  7c  inner  scale  [m] 

T .  Lo  =  inf;  7,  outer  scale  [m] 

T. Model  =  ’FHill’;  7c  spectral  model 

T.fmin  =  min(Fx); 

T.fmax  =  max(Fx); 

T.N  =  2  ~7 ; 

T  .  rO  =  rO  ; 
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7c  min  spatial  frequency  [1/m] 

7c  max  spatial  frequency  [1/m] 

7c  number  of  FS  coefficients 
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7.  loop  over  realizations 

phi  =  zeros (N  ,  N)  ; 

for  idx.num  =  1  :  numscreens 

idxStr_num  =  spr intf  (  1  7«03d  ’  ,  idx.num)  ; 


PSD  =  mkPSD_thesis  (T)  ;  7,  create  the  PSD 

C  =  FSCoef  f  _the  s  i  s  (  PSD  )  ;  7,  generate  coefficients 

x  =  (-N/2  :  N/2-1)  *  deltax;  7«  position  coordinates  [m] 
y  =  x; 

S  =  FSscrnREAL_thesis  (x  ,  y,  C);  7«  generate  the  phase  ... 
screens 

phi ( : , : )  =  S . g-mean (mean (S . g) ) ; 


l  Calculate  zernike  coefficients  to  for  Piston,  Tip,  and  . 
Tilt 

a.piston  =  sum ( sum (circ_ap.*phi(:  ,  :)  ,*Z_piston))  ./ sum ( 
sum(circ_ap ,*Z_piston  .  ~2) )  ; 

a_xtilt  =  sum ( sum (circ_ap.*phi(:  ,  :)  ,*Z_xtilt))  . / sum ( . . 
sum(circ_ap ,*Z_xtilt  .  ~2) )  ; 

a_ytilt  =  sum ( sum (circ_ap.*phi(:  ,  :)  ,*Z_ytilt))  . / sum ( .  . 
sum(circ_ap ,*Z_ytilt  .  ~2)  )  ; 

70  Remove  Piston,  Tip,  and  Tilt  from  the  phase  screens 

phi_PTR(:,:)  =  phi ( : , : )  -  (a_piston.*Z_piston)  -  (... 

a_xt ilt . * Z_xt ilt )  -  ( a_yt ilt . * Z_yt i It )  ; 


7»  Save  important  variables  for  future  use 

save ( [ ’Thesis\Matlab\mat  Files\thesis_DrO 
’ .screens. ’  idxStr.num  ’  . mat  1 ]  ,  .  .  . 

’ phi ’  ,  1  phi _PTR ’  ,  ’D’,  ’D.over.rO’  ,  ’ 

,  'theta’,  ’numscreens'); 


idxStr _DrO 


circ.ap  ’  , 


end 


Appendix  C.  Matlab  Code  to  Pixelate  and  Compute  Statistics  of 

Phase  Screens 

Listing  C.l:  Matlab®  file  to  down-sample  and  compute  statistics  of  phase  screen 
realizations. 


7,  Fourier  Series  Phase  Screen  Evaluator 
clear  all ;  clc  ; 
format  compact; 


5  7«  load  path  for  locations  of  functions  called 

addpath (’ Thes is \ Mat  lab \ Final  Thesis  Files’); 
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15 


7,  Define  important  parameters 

numscreens  =  250;  7«  number  of  random  phase  screens  to  create 

N  =  512;  7.  number  of  screen  grid  points 

D  =  2;  7.  aperture  diameter  [m] 

D_over_r0  =  5  7«  strength  of  turbulence  to  be  created 

idxStr_DrO  =  sprintf  (  1  7«02d  ’  ,  D_over_r0); 
cor_Z_mode_PLOURDE_percent_pres_mode_equation  =  139;  °i  ... 
number  of  modes  required  for  threshold  level 

7.  Create  array  of  values  for  N_new 

pixels  =  [512  256  128  64  32  16  8  4  2  1] ; 


7.  Create  a  circular  aperture  of  unit  radius 
20  x  =  linspace (-D/2 , D/2  ,  N); 

[x  y]  =  meshgrid (x , x)  ; 

[theta  rho]  =  cart2pol (x , y ) ; 
circ_ap  =  rho  <  D/2; 
r  =  (2 . *  rho )  . /D ; 

25 

7.  Pre-compute  all  of  the  Zernike  polynomials  using  eval.zernike  ... 
function 

for  idx_mode  =4  :  ... 

cor_Z_mode_PLOURDE_percent_pres_mode_equation 
Z_mode(:,:,idx_mode)  =  eval_zernike(idx_mode,  r  ,  theta); 
denom(idx_mode)  =  sum( sum (circ_ap.*Z_mode(:  ,  idx.mode) .  .  . 

•  ~2))  ; 

30  Z.mode , idx.mode)  =  circ_ap.*Z_mode(:,:,idx_mode); 

end 


7.  Iterate  over  all  pixel  at  ion  levels 
for  idx_pixels  =  1 : length (pixels ) 

35  idxStr_pixels  =  spr intf  (  ’  7o03d  ’  ,  pixels  (  idx_pixels  ))  ; 

iter  =  pixels ( idx_pixels ) ; 

idxStr_iter  =  sprintf  (’ 7«03d  ’  ,  iter) 

7,  Iterate  over  all  random  phase  screens 
40  for  idx.screen  =  1  :  numscreens 

idxStr.screen  =  sprintf  (’ 7003d  ’  ,  idx_screen); 


90 


45 

50 

55 

60 

65 

70 

75 


7.  Load  each  PTR  phase  screen  realization 

load ( [ ’ Thesis\Matlab\mat  Files\thesis_DrO_ 
idxStr.screen  ’ .mat ’  ] , . . . 


idxStr _DrO , 


.screens. 


’ phi_PTR 1 ) ; 

if  pixels ( idx_pixels )  ~=  512; 

/  Down-sample  from  N  to  N/div 

pix  =  pixels  ( idx_pixels )  70  how  many  pixels  you 

want  to  end  up  in  each  row 
N  =  length  (phi_PTR)  ;  70  how  many  pixels  you 

are  starting  with  in  each  row 
div  =  N/pix  ;  70  the  reduction  in  ... 

pixels 


for  idx  =  l:N/div; 
pull  from 
odd(idx)  =  count; 
count  =  count  +  div; 

end 


7,  create  a  matrix  to 


for 


end 


idx2  =  1 : N/div 
for  idx  =  l:N/div 

odd_val ( idx2 , idx )  =  mean ( mean ( phi_PTR ( odd ( . 
idx2)  :odd(idx2)+(div-l)  ,odd(idx)  : odd ( . .  . 
idx  )  +  (  div  -  1 )  )  )  )  ; 


1 


Create  pixelated  PTR  phas 
realization 


screen 


phi_PTR_pixellated(odd(idx2) : odd(idx2) 
+  (div-l)  ,  odd(idx)  : odd(idx)  +  (div-l) 
)  =  odd_val ( idx2 , idx )  ; 


end 


end 


’/,  Calculate  Zernike  coefficients  for  the  pixelated  PTR  ... 

phase  screen  realization 
for  idx.mode  =  4  :  ... 

cor_Z_mode_PLOURDE_per cent _pr es_mode_equat ion 
if  pixels ( idx_pixels )  ~=  512; 

a_mode ( idx_mode ,  idx.screen)  =  sum(sum(... 
phi_PTR_pixellated(:  ,:)  . *Z_mode(:  ,  ... 

idx_mode))) ,/denom(idx_mode) ; 

else 

a_mode ( idx_mode ,  idx.screen)  =  sum(sum(... 
phi_PTR(: , :) ,*Z_mode(: , : ,  idx_mode))) ./. 
denom ( idx_mode  )  ; 

end 

end 


end 


7„  Comput e  the  statistics  of  the  PTR  phase  screen  realizations 
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i  Find  variance  in  each  Zernike  mode 
phvar_a_mode  =  ?ar(a_mode  ,0,2)  ; 

80 

7.  Find  HO  total  phase  variance  up  to  each  Zernike  mode 
HO_cum_phvar_a_mode  =  cumsum ( phvar_a_mode ) ; 

%  Convert  to  phase  standard  deviation 
85  HO_cum_phstdev_a_mode  =  sqrt (HO_cum_phvar_a_mode ) 

.  / (2*pi)  ; 

it  Remove  r_0  dependence  from  the  phase  standard  deviation 
HO_norm_cum_phstdev_a_mode  =  HO_cum_phstdev_a_mode  . / . 
D_over_r0  .  “(5/6)  ; 

90  %  Find  total  HO  phase  std  dev  up  to  each  Zernike  mode  for 

each  pixelation  level 

HO_norm_cum_phstdev_a_mode_ALL ( idx_pixels +  1  ,  : )  =  ... 

HO_norm_cum_phstdev_a_mode ; 

‘/,  Find  maximum  possible  phase  standard  deviation  for  ... 
simulation  and  show  on  main  screen 
phstdev_pres_af ter_modes_and_pixel  =  ... 
HO_norm_cum_phstdev_a_mode  (  . .  . 

cor _Z_mode_PLOURDE_per cent _pres_mode_equat ion ) 

95  total_phstdev_pres_af ter_pixel  =  max ( . .  . 

HO_norm_cum_phstdev_a_mode ) 

*/,  Create  array  of  pixelation  levels  and  the  associated  .. 
simulation  phase  standard  deviation  values 
Computer_SIM_phstdev_pres_after_modes_and_pixel  (. .  . 
idx_pixels+2 , 1)  =  iter; 

Computer_SIM_phstdev_pres_after_modes_and_pixel  (. . . 
idx_pixels +2 , 2)  =  ... 

phstdev_pres_after_modes_and_pixel; 

100 

*/,  Create  sample  PTR  phase  screen  realizations  for  figures 
phi _PTR_ samp (:,:,idx_pixels+l)  =  phi_PTR; 
phi _PTR_pixe 1 1  at ed_ samp  (  :  ,  :  , idx_pixel s +1 )  =  ... 

phi_PTR_pixellated ; 

105  7.  Save  important  variables  for  future  use 

save ( ’ Thesis\Matlab\mat  Files \H0_computer_SIM_re suits . mat 

J  •  •  • 

’ Computer_SIM_phstdev_pres_after_modes_and_pixel  ’  ,  ’  . 

phi_PTR_samp ’  ,  .  .  . 

’ phi_PTR_pixellated_samp  1  , 

HO_norm_cum_phstdev_a_mode_ALL 1 ) 

end 
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